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ABSTRACT 


Pearce’s  (1977)  processor  for  analyzing  seismic  focal  mechanisms  has  been 
implemented  on  two  seismic  analysis  systems,  and  it  has  been  given  limited 
testing  for  use  as  a  discrimination  aid.  This  algorithm  uses  relative  amplitudes 
of  the  pP  and  sP  depth  phases  with  respect  to  the  P- wave  amplitude  at  each 
station  in  a  network  If  pP  and/or  sP  cannot  be  identified,  then  the  amplitude 
of  the  P-wave  coda  at  the  anticipated  arrival  time  of  the  undetected  depth 
phase  may  be  used  as  an  upper  limit  on  the  amplitude  of  that  phase.  The  solu¬ 
tion  space  of  double-couple  seismic  sources  is  searched  for  possible  solutions 
consistent  with  the  relative  amplitude  bounds  which  were  measured  at  each  sta¬ 
tion.  If  no  solution  can  be  found,  then  either  the  "depth  phases”  were 
misidentified  or  the  event  was  not  a  double-couple  (i.e. ,  it  may  have  been  an 
explosion)  The  algorithm  and  support  software  have  been  installed  at  the 
Seismic  Research  Center  (SRC)  on  the  Seismic  Research  Information  System 
(SRIS)  and  at  the  Center  for  Seismic  Studies  (CSS).  A  description  of  both  sys¬ 
tems  is  given.  Examples  of  the  use  of  the  algorithm  as  a  seismic  analysis  tool 
are  presented. 

Theoretical  analysis  has  been  performed  to  evaluate  the  utility  of  the  con¬ 
straint  on  the  amplitude  of  the  depth  phases  which  is  provided  by  the  amplitude 
of  the  seismic  coda  as  a  discrimination  tool  for  testing  the  hypothesis  that  the 
event  is  an  earthquake  at  depth.  Analysis  using  synthetic  data  suggests  that  the 
Pearce  algorithm  may  be  a  useful  discriminant  of  explosions  when  the  coda 
remains  below  the  initial  P- wave  amplitude  at  five  or  more  teleseis'mic  stations 
covering  all  four  quadrants.  In  such  a  case,  it  is  highly  probable  that  the  algo- 
rithm  will  not  find  any  double-couple  solutions  at  depth  consistent  with  small  pP 
amplitudes  at  all  stations.  A  test  on  Global  Digital  Seismic  Network  (GDSN)  data 
suggests  that  pP  amplitude  constraints  are  sufficient  to  eliminate  the  possibility 


of  double-couple  sources  at  depths  greater  than  3-5  km  for  explosions  of  m*  >  5, 
when  azimuthal  coverage  is  adequate.  Testing  of  the  proposed  discriminant  on 
events  with  smaller  magnitudes  could  not  be  adequately  performed  with  GDSN 
data.  Analysis  of  Eurasian  earthquake  GDSN  data  demonstrated  thatpP  phases 
could  be  detected  on  short  period  records  for  shallow  (less  than  33  km)  events, 
with  acceptable  double-couple  solution?. 
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L T ST  OF  FTGFRFS 


T  i  1 1  e  Pag 

Figure  6  reproduced  from  Pearce  (1979),  definitions  of  1 

the  source  parameters  strike,  slip,  and  dip  for  the 
fault  plane,  and  the  takeoff  angles  for  P,  pP ,  and  sD. 

Figure  7  reproduced  from  Pearce  1197*1,  showing  the  1 

representation  of  the  three  dimensional  closed  solution 
space  of  doub  1  e- c ou p 1 e s  represented  hv  slip,  dip,  and 
strike  angles  used  in  the  Peaice  algorithm. 

P  wave  amplitude  and  coda  levels  used  as  input  to  the  0 

focal  plane  algorithm  as  p,  pP ,  and  sP  amplitude  hounds 
are  indicated  with  horizontal  bars  (event  149  Table  T.) 

F  vent  149  10/12/80  ■  9  :  24  : 1 4  .  1  0  ,  49.96':,  7O.08F.  0  km,  ? 

m^5.9  M  4.1.  (h)  A  histogram  of  the  percentage  of  accept¬ 
able  doub 1 e- c oup 1 e s  is  shown  for  the  number  of  incon¬ 
sistent  observations  tolerated.  \  total  of  20  observations 
for  the  10  stations;  pP/P  and  sP/P  bounds  for  each  station. 
Oumo  and  0 r f o  were  not  used.  Three  of  the  observations 
must  he  discarded  before  solutions  are  allowed.  (c)  An 
equal  area  lower  hemisphere  projection  of  the  p  wave  take 
off  angles  for  each  station.  Tf  Kono,  Anto,  and  Bcao  are 
discarded  then  two  quadrants  are  open  and  solutions  are 
admitted  by  the  focal  piane  algorithm. 

Seismograms  for  the  'lew  Brunswick  earthquake  adapted  2 

from  Chov  et  al  (1981).  The  estimated  P,  pD  and  sP 
amplitude  hounds  are  indicated  on  the  short  period 
seismograms . 

Acceptable  focal  mechanisms  for  the  'lew  Brunswick 
earthquake.  Solutions  span  a  space  of  d ou b 1 e- c ou p 1 e s 
with  well  constrained  tension  axis  and  a  variable 
compression  axis.  Tension  (")  and  Compression  (P)  axis 
are  shown  on  a  equal  area  lower  hemisphere  focal  sphere. 

P  wave  take-off  angles  are  plotted  for  each  station  with 
the  station  name.  P  wave  nodal  planes  are  shown  for  each 
acceptable  solution. 

Several  solutions  from  Chov  et  al  (1°81)  and  Hzievonski 
and  Woodhouse  (1981). 

Seismograms  reproduced  f  r  m  Pearce  (1979)  for  the  May  1, 
1969  event. 

\cceptable  solutions  to  t  »>e  amplitude  constraints  shown 
in  (A)  . 
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Title  v 


Acceptable  solutions  for  the  March  20,  1976  event  using 

only  pP/P  and  sP/P  amplitude  ratio  constraints.  Fach  set 
of  P  and  T  axis  is  exchanged  because  polarites  were  not 
considered. 

Simple  one  layer  crustal  model  synthetics  for  three 
different  possible  solutions. 

Simple  one  layer  crustal  model  synthetics  for  three 
different  possible  solutions. 

Simple  one  layer  crustal  model  synthetics  for  three 
different  possible  solutions. 

I'nfiltered  seismograms  of  event  1716. 

T'nfiltered  seismograms  of  event  1914. 

Filtered  seismograms  of  event  1914. 

b'FTS  reported  polarities  for  1716. 
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algorithm. 
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planes 
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from 

the 

Pearce 

A  r  r  e  p  t  a  b  1  c 
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1914 

from 

"nf  i  I  1 1  ri'd  seismograms  of  event  2f>. 

Filtered  seismograms  of  event  20. 

1  ■  ■  reported  polarities,  and  acceptable  solutions  iron 
the  Pearce  algorithm. 

Seismograms  for  event  R  on  the  Siberian  Coast. 

Acceptable  double-couples  for  event  22  consistent  with 
pP/P  and  s P/P  bounds  at  Ion,  jas,  alq,  bcao,  ctao,  and 


ETST  ')!•'  F  T  fl  I TR  F.  S  (Continupd) 

Figure  No.  Title  Pape 

He  Polarity  data  for  event  ?2.  Note  the  cluster  of  6  4 

compressions  in  the  NW  quadrant  apparently  inconsistent 
with  the  solutions  in  8b. 

9a  Seismograms  for  event  8  on  the  Siberian  Coast.  5 A 

9h  Acceptable  double-couples  for  a  set  of  tight  bounds  on  5 7 

sP,  on  the  left.  4  observations  have  been  ignored  to 
generate  this  solution  set.  TSC  reported  polarities  are 
plotted  on  the  right. 

9c  Double-couple  solution  set  for  a  sc^t  of  relaxed  sp  58 

amplitude  hounds .  On  the  left,  P,T,  and  B  axis  without 
the  clutter  of  focal  planes.  On  the  right,  the  clutter 
of  focal  planes. 

If  a  Seismograms  For  event  295.  Predicted  arrival  times  are  An 

shown  for  the  P,  PcP,  pP,  and  sP  times  for  the  NFTS 
hypocentral  location. 

10b  Synthetics  using  the  simple  crustal  model.  The  hi 

algorithm  does  not  properly  show  depth  phases  for  the 
Pk?  arrivals  at  Anmo  and  Zobo.  Amplitudes  are 
normalized  to  maximum  possible  P  amplitude  from  a 
double-couple  source  at  that  distance.  Strike,  slip  and 
dip  from  NF.TS.  Note  the  large  pP  and  s?  arrivals  with 
respect  to  the  P  wave  amplitude. 

10c  NF.TS  double-couple  solution  with  reported  NFTS  6? 

polarit  ies. 

lOd  Acceptable  solutions  from  the  Pearce  algorithm  based  on  68 

pP/P  amplitude  ratio  bounds  only.  Two  subsets  of  the 
solution  space  are  apparent. 


lOe 

!  1  a 

i  lb 

1  1  c 


Acceptable  solutions  from  the  Pearce  algorithm  based  on  f>  A  -i  j 

pP/P  amplitude  ratio  bounds  only.  Two  subsets  of  the  jj 

solution  space  are  apparent.  ! 

Seismograms  for  earthquake  516.  65 

NETS  polarities  for  earthquake  516,  with  the  NFTS  6n 

doub 1 o - c ou p 1 o  .  Note  the  large  number  of  inconsistent 

polarities.  | 

Acceptable  solutions  from  the  Pearce  algorithm  analysis  67 
of  the  data  in  (A).  Polarities  were  onlv  used  for  Col, 

To  1  ,  and  Nwao.  Bcao,  Humo,  Ctao  and  Zobo  appear  to  show  I 

nodal  P  waves  In  (A).  i 
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F igure  No . 


Title 


12a 


12b 


12c 


11a 


lib 


14 


The  percentage  of  all  possible  double-couples  eliminated 
by  the  constraint  of  |pP/Pj  <  1 ,  and  |sP/P|  <1  at  a 
single  station.  Constraints  are  considered  together  and 
separately.  The  teleseismic  takeoff  angle  is  varied 
from  10  to  25  degrees.  The  observation  that  |pP/P|  <  1 
eliminates  nearly  50?  of  all  double-couples  at  any 
teleseismic  takeoff  angle.  The  implications  that  sP  P 
are  less  significant  and  slightly  dependent  on  takeoff 
angle.  The  observation  that  both  pP  and  sP  are  smaller 
than  P  implies  that  slightly  fewer  solutions  are 
acceptable  the  observation  the  pP <  P  alone. 

The  percentage  of  all  possible  double-couples  eliminate 
by  one,  two,  three,  or  four  stations  at  10  degree 
azimuthal  increments  varies  as  a  function  of  the  absol 
limit  of  the  pP  and  sP  amplitudes.  All  takeoff  angle? 
are  assumed  to  he  20  degrees. 

The  percentage  of  all  possible  double-couples  eliminate 
by  a  network  of  stations  at  60  degree  azimuthal 
increments.  Four  or  more  stations  are  required  with 
|pP/p|  and  |sP/P|  <  1  to  eliminate  all  possible  focal 
mechanisms  from  consideration. 

Results  of  a  hypothetical  nine  station  network  with  two 
arrivals  of  nearly  equal  amplitude.  Both  P  and  the 
prospective  pP  are  given  25?  amplitude  tolerances.  462< 
possible  focal  planes  out  of  12212  are  found.  No 
polarity  information  is  used.  Solutions  are  either  of 
dip-slip  with  normal  thrust,  dip-slip  with  reverse 
thrust,  vertical  dip-slip,  or  horizontal  faults  with 
horizontal  slip. 

Results  of  a  hypothetical  nine  station  network  with  two 
arrivals  of  nearly  equal  amplitude.  Both  P  and  the 
prospective  pP  are  given  257,  amplitude  tolerances.  462 
possible  focal  planes  out  of  12112  are  found.  No 
polarity  information  is  used.  Solutions  are  either  of 
dip-slip  with  normal  thrust,  dip-slip  with  reverse 
thrust,  vertical  dip-slip,  or  horizontal  faults  with 
horizontal  slip. 

Flow  chart  of  data  analysis. 


15  P-wave  amplitude  bounds  measured  using  the  AMPT.TM 

opt  ions . 
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Table 
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TT 

ITT 


Mo.  Title 

Explosions  in  the  OSS  data  base  used  for  testing  the 
Pearce  focal  plane  analysis  algorithm.  Depth  fh) 
constrained  to  zero. 

Earthquakes  added  to  the  OSS  data  base  for  testing  of 
the  Pearce  focal  plane  analysis  algorithm.  Depth  ih) 
in  km  unconstrained  unless  h  =  'll  km. 

SIJBMKNI’  Options  of  the  pearce  Measurement  Subsystem 


TV 


Polynomial  pit  to  the  Tables  of  P-wave  "rake-off  Angles 


INTRODUCTION 

The  traditional  method  of  determining  the  focal  mechanism  cf  ari  earth- 
quake  has  been  to  plot  on  ,»  stercographie  projection  the  polarities  of  !  he  first 
motions  of  the  P  waves  which  were  observed  for  that  earthquake  at  a  net  ac  rk  of 
stations  An  orthogonal  pair  of  large-circle  arcs  could  then  be  found  graphically 
which  best  separated  the  eomprcssicnal  and  dilational  first  motion-  into  qua¬ 
drants.  These  arcs  represent  the  projection  of  the  fault  and  auxiliary  planes  of 
the  double-couple  mechanism,  and  their  strike  and  dip  could  be  read  directly  ofT 
the  stereonet  The  rake  of  the  slip  vector  could  also  be  calculated  by  not.ng 
that  the  slip  vector  is  perpendic  ular  to  the  intersection  of  the  fault  and  auxili¬ 
ary  planes  (i  e  ,  the  null  axis),  so  all  the  parameters  necessary  for  specifying 
the  focal  mechanism  could  be  determined  from  this  graphical  srpar  d-rn  of  the 
first-motion  polarities  into  quadrants  Kvcntu.illy  computer  searches  ever  all 
possible  focal  mechanisms  were  used  to  replace  the  sometimes  subject ivo  pro¬ 
cess  of  physically  manipulating  a  stereonet  in  order  to  find  the  best-filling  pair 
oi  orthogonal  planes  (e  y  ,  Wickens  and  Hodgson,  ‘.967) 

Although  this  traditional  technique  is  sufficient  in  principle  tc  permit  the 
determination  of  the  earthquake  focal  mechanism,  it  uses  only  a  fraction  of  the 
information  pertaining  to  the  seismic  radiation  pattern  which  is  available  for 
analysis  from  the  seismogram  The  traditional  technique  has  therefore  been 
expanded  in  the  last  20  years  or  so  to  incorporate  additional  information  other 
than  simply  the  P- wave  polarities.  One  such  type  of  additional  information  is 
the  polarity  of  the  5- waves,  which  can  be  measured  by  taking  the  ratio  of  the 
amplitudes  of  the  Slf  and  .ST  polarizations  The  5 -wave  radiation  p.Hc'o  is 
such  that,  if  the  directions  of  the  .^-wave  polarities  arc  plotted  as  arrows  on  the 
same  stereographic  projection  as  was  constructed  Tor  the  P- wave  polarities,  the 
arrows  will  be  directed  from  the  compression  axes  to  the  tension  axes,  and  they 
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ni  after  the  ••winj*  of  the  first  motion,  then  rort-un  a  number  of 
tr-  this  augmented  approach  !-’rr  one  thing  tv  ...r-sourte  focusing 
»!  irrionuigeneitirs  along  the  prrp.14.it ion  path  may  cause  il'.c  sect  >n 
•nine  radiation  pattern  which  is  actually  sampled  by  a  giv-'n  station  to 
n:  : r  r  m  iluit  which  would  bo  predicted  from  the  a/ir  adh  and  take-off 
as  :,r  1.  d  m  plot ,ng  lh..i  status.  c  n  in«  mb-'  ma\ 

a  t  i-ither  on  account  of  poor  signni-t o-noisr  r.-.t'o  or  o-  a- a  c  i"t  of  a 
incri-.o-inj;  seismic  source  time  function,  and  these  emi  >s.  -  i  signab 


cnlogori/ed  on  the  basis  of  their  polaril j  Furthermore  the  choice  of 


the  particular  time  window  in  which  one  should  measure  the  amplitude  of  an 
emergent  signal  may  be  rendered  ambiguous.  Focusing  near  the  receiver  and 
uncertainty  in  the  value  which  is  assumed  for  the  anelastic  attenuation  parame¬ 
ter  t  *  for  the  souree-to-receiver  propagation  path  may  result  in  the  over-  or 
underestimation  of  the  signal  amplitude  of  either  type  of  body  wave  Calculation 
of  the  seismic  moment  tensor  will,  of  course,  be  particularly  sensitive'  to  errors 
in  the  amplitude  measurement  or  the  path  correction  for  stations  close  to  the 
nodal  planes  of  the  double  couple  Finally,  5-waves  cannot  be  expected  to  con¬ 
tribute  much  information  to  the  solution  for  small  events 

Pearce  (1977,'. 9t30)  has  developed  an  e  xtension  to  the  traditional  method  of 
focal  mechanism  solution  which  offers  certain  advantages  over  the  methods 
described  above  Pearce’s  technique  relies  upon  the  formation  of  the  ratios  of 
the  observed  amplitudes  of  the  phases  P,  pP.  and  sP  Telcscismic  radiation 
from  the  surface-reflected  phase  p!}  samples  the  upper  half  of  the  focal  sphere 
which  is  not  sampled  by  the  downward  propagating  P  waves,  so  the  chance?  of  a 
given  station’s  being  close  to  a  nodal  plane,  and  hence  of  exerting  strong 
influence  on  the  fault  plane  solution,  are  increased  when  pP  amplitudes  arc 
added  to  the  data  base  of  P  amplitudes  Use  of  the  sP  phase  allows  one  to  sam¬ 
ple  the  5-wave  radiation  pattern  without  having  to  measure  5  directly  Perhaps 
most  importantly,  forming  the  ratios  P/pP,  P/sP,  and  pP/sP  tends  to  cancel 
out  the  effects  of  attenuation,  since  the  propagation  paths  for  each  of  these 
phases  are  nearly  the  same  for  shallow  events.  By  using  sP  rather  than  the 
direct  5  phase,  one  can  ignore  the  effects  of  t  * s  on  all  of  the  path  except  for  the 
upgoing  leg,  so  the  amplitude  ratio  /V sP  which  is  used  by  Pearce  is  less  sensi¬ 
tive  to  the  effects  of  attenuation  than  is  the  ratio  P /  5,  which  was  used  for  focal 
mechanism  solutions  by  Honda  (1957)  Another  important  advantage  of 
Pearce's  technique  is  that,  if  pP  or  sP  cannot  be  detected,  the  amplitude  of  the 
seismic  coda  at  the  time  of  the  predicted  phase  arrival  can  be  measured,  and 
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this  value  can  then  be  used 


u  e  i.ar'a')- 


<is  an  upper  hr  and  ' 

served  phase. 

In  order  to  implement  rearer  wv.pht  urn-  r.t* 
to  have  an  earth  model  fee  ray  tracing.  which  ;>«  r::.r.  -  1  '  ■  •:  -  •'  •  ate  '  r.e 

measurements  made  at  a  pven  station  w  '  h  tin  :••••;  :••<  ;  •  • '  n  <  :  : • 

sphere  from  which  the  seismic  radiation  erninatc.-a  u  *v  r<  •  .  • :  <_•  same 

as  the  traditional  technique  which  as-o  rates  st  at  ions  t.:,  ‘lie  st<  rr  np.i  t  with 
points  on  the  /'-wave  radiation  pattern  ov  means  of  • ...  a.'.m.uth  and  take-off 
angle  One  further  piece  of  information  :>  required  fur  Pearce  s  metr.ud.  how¬ 
ever.  namely  an  eart  h  mode!  which  accounts  for  the  npgoing  and  r  '-fleeted  leg.-: 
(back  down  to  the  soun  e  dept h)  of  pP  and  s/’  !n  the  c  asc  of  pP.  a  simplifying 
issumption  is  made  namely  that  the  event  is  at  shallow  depth  so  the  only  effect 
which  the  la \ ers  nf  the  earth  mode!  above  the  source  depth  have  on  the  propa- 
gat  .on  of  />/•'  is  to  introduce  a  simple  multiplicative  cocfTu  lent  representing  the 
efficiency  of  the  reflection  at  the  free  surface  in  the  case  of  <?/•*.  it  is  also 
assumed  that  the  event  is  shallow,  and  the  effect  of  ■  hr  shallow  layers  is  to 
refract  the  upgning  .‘'’-wave  leg  and  reflected  /’-wave  leg  in  a  simple  manner 
which  can  be  calculated  using  only  the  seismic  velocities  at  the  source  depth 
and  at  the  free  surface,  the  algorithm  thus  uses,  in  effect  a  one-layer  crustai 
model  Again,  a  coefficient  is  also  introduced  which  represents  I  he  efficiency  of 
the  surface  reflection,  this  time  taking  into  account  the  conversion  from  to 
/'-wave  energy  We  shall  discuss  the  implications  of  thest  amplifying  assump¬ 
tions  and  the  choice  of  the  reflection  coefficients  later 

It  can  be  seen  immediately  that  one  difficulty  with  I ’carte  s  technique  is  t  h» 
use  of  this  simplistic  model  for  the  propagation  and  i  cfVt  •  a  ■  ,  I  -  c  :;v.o:u_ 
legs  of  pP  arid  sP  There  are  also  certain  practical  problems  >n  m  ikmg  the 
necessary  measurements  of  the  sought-after  phases  on  Inc  sciMin.gram.  a.-  we 
shall  point  out  later  It  is  also  important  lo  note  that  Pearce's  technique  doe  - 


not  escape  certain  of  the  difficulties  which  beset  the  other  extensions  of  the 
traditional  fault  plane  solution  which  were  discussed  previously.  Specifically, 
near-source  focusing  and  lateral  inhomogeneities  along  the  propagation  path 
can  cause  the  portion  of  the  radiation  pattern  sampled  at  a  given  station  to  be 
different  from  that  which  is  expected  on  the  basis  of  the  azimuth  and  take-off 
angle  which  were  used  to  plot  that  station  on  the  stereonet;  this  problem  is  now 
compounded  by  similar  problems  for  the  reflected  legs.  Near-receiver  focusing 
will  still  be  a  source  of  error,  since  this  may  fail  to  cancel  out  when  the  ampli¬ 
tude  ratios  are  formed  of  phases  which  have  reached  the  receiver  region  by  even 
slightly  different  paths.  Emergent  signals  will  still  contribute  no  useful  polarity 
information,  and  it  may  be  difficult  to  determine  where  in  the  wavetrains  of 
these  signals  the  measurements  of  P  and  of  the  bounds  for  pP  and  sP  should  be 
made. 

We  have  mentioned  that  amplitude  information  for  the  P  and  S  phases  is 
introduced  into  the  focal  mechanism  solution  by  means  of  inversion  to  deter¬ 
mine  the  seismic  moment  tensor.  By  contrast,  in  the  implementation  of 
Pearce's  technique  amplitude  information  is  introduced  by  means  of  a  grid 
search  over  all  possible  geometrical  orientations  of  the  double  couple.  These 
orientations  are  defined  by  allowing  the  strike  of  the  fault,  the  dip  of  the  fault, 
and  the  rake  of  the  slip  vector  within  the  fault  to  take  on  all  possible  values.  In 
distinction  to  the  more  usual  convention  of  defining  the  double-couple  orienta¬ 
tions  by  0°Sstrike<360° ,  0®sdip<900,  and  — 180°^rake<180° ,  Pearce  defines 
them  by  0,^strike<360o,  0°;sdip<1800 ,  and  0°^rake<lB0° .  The  definitions  of 
each  of  these  parameters,  including  the  sign  conventions,  should  be  inferred 
from  Figure  1.  We  shall  perform  the  grid  search  over  strike,  dip,  and  rake  by 
sweeping  over  each  of  these  three  parameters,  within  the  limits  listed  above,  in 
increments  of  either  5#  or  10*.  In  presenting  graphically  the  results  of  his 
analysis,  Pearce  plots  each  allowed  focal  mechanism  as  an  individual  vector 


ominaling  from  the  corresponding  point  within  the  plane  spanned  by  the  coordi¬ 
nates  dip  and  rake,  and  having  as  the  vector's  azimuth  the  strike  »  '  tut  double 
couple  (c/  Figure  2). 

For  every  point  on  the  strike-dip-rake  grid  defined  by  to.-  ! . , t  -■  ’ist.ed 
above,  the  ratios  P/pP,  P/sP,  and  pP/ sP  are  computed  "or  '  voy  station  c 
the  network.  These  values  are  then  compared  with  the  ratios  of  the  observed 
amplitudes  and  with  the  observed  polarities.  (If  no  amplitude  >r  format, on  ever, 
upper  bounds  deduced  from  the  coda  amplitude,  is  available  for  a  giver,  phase  at. 
o  given  station,  then  the  ratios  involving  tnat  phase  ae  simply  i gne-ed  This  also 
rppUos  to  the  polarity  of  each  of  the  three  phases  )  Of  course-  the  predicted 
amplitude  ratio  should  not  be  expected  to  be  exactly  equal  !.o  the  observed 
value,  so  it  is  necessary  to  specify  upper  and  lower  bounds  abound  the  observed 
■-mphtudes  of  each  phase  representing  itv  error  bars"  v.hc  h  t*  •n,r.  -u  ,-<j 
geophysical  effect  s  such  as  focusing  as  well  as  the  imprecision  of  invasur  truml 
In  tne  ease  of  an  unobserved  pP  or  sP  phase,  the  lower  bounu  u>  taaen  ic  be 
zero,  and  the  upper  bound  is  taken  to  be  the  observed  coda  amplit  ude  or  a 
slightly  larger  value  than  that  Instead  of  forming  the  ratios  of  the  observed 
amplitudes,  then,  Pearce's  algorithm  actually  forms  ratios  cf  the  amplitude 
bounds i,  and  these  are  in  turn  used  as  bounds  on  the  observed  ratio:-.  What  is 
tested  is  thus  whether  the  predicted  amplitude  ratios  fall  within  the  computed 
bounds  around  the  observed  ratios,  and  whether  the  predicted  polarities  agree 
with  those  which  are  observed.  If  this  comparison  fails  for  at  least  one  station  in 
the  network,  then  this  grid  point  is  eliminated  from  the  suite  of  possible  focal 
mechanism  solutions. 

The  grid  search  actually  proceeds  by  testing  the  predicted  amplitude  ratios 
and  polarities  for  each  grid  point  against  the  observed  bounded  ratio  and  polart- 
t'es  of  a  given  phase  pair  at  a  given  station,  and  then  the  search  over  t  he  entire 
grid  is  repeated  for  the  next  phase  pair  at  that  station  or  for  a  phase  pair  at  the 
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Figure  7  Method  of  representing  acceptable  fault  plane 
orientations  in  terms  of  slip  direction  o ,  dip  -s  and  strike  as 
defined  in  Figure  6(a).  Acceptable  orientations  are  plotted  as 
vectors  from  the  Cartesian  point  defining  u  and  <'>,  in  the 
direction  of  the  strike  u.  Lower  hemisphere  stereographic 
projections  indicate  the  type  of  fault  plane  orientation  represented 
by  various  combinations  of  and  ,1,  and  are  shown  oriented  for 
strike  =  360°  (northerly) 

Figure  2  Figure  7  reproduced  from  Pearce  (1979,  showing  the 

representation  of  the  three  dimensional  closed  solution  space 
of  double-couples  represented  by  slip,  dip,  and  strike  angles 
used  in  the  Pearce  algorithm. 
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nr  xt  stal. on  i  r ;  1 :  ill  t  he  observed  data  have  been  tested.  A  list  of  the  grid 
noints  (•-  a  the  double  couple  mechanisms)  which  are  consistent  with  ihe  par¬ 
ticular  riati-.  point  under  consideration  is  printed  out  along  with  a  cumulative  list 
of  the  grid  points  which  are  consistent  with  all  data  tested  so  far.  At  i.he  end  of 
the  process,  t.i.s  cumulative  list  contains  all  of  the  double-couple  mechanisms 
which  are  considered  to  be  acceptable  solutions.  Note  that  this  list,  will  almost 
always  encompass  a  range  of  solutions  rather  than  a  unique  value,  such  as  is 
usually  presented  to  be  "the"  fault  plane  soltuion  which  is  found  using  polarities 
plotted  on  a  stereonet  Note  also  that  the  range  of  solutions  cannot  always  be 
represented  as  a  single  value  with  error  bars  around  each  of  tne  three  double- 
couple  parameters,  since  the  grid  points  which  satisfy  ail  the  tests  are  not 
necessarily  contiguous.  An  extension  to  Pearce's  algorithm  which  we  have 
implemented  in  this  study  is  to  generalize  this  logic  to  pet-mi  the  'el.nt  on  of 
grid  points  for  which  the  number  of  inconsistent  observations  is  ie^s  than  some 
specified  limit,  rather  than  automatically  rejecting  them  if  the  nurnbet  of  incon¬ 
sistencies  exceeds  zero  This  revision  prevents  the  solution  from  being  dom¬ 
inated  by  a  small  number  of  possibly  erroneous  measurements.  (Incidentally, 
we  note  that  one  common  source  of  measurement  error,  namely  miscall  oration 
of  the  instrument  gain,  is  inapplicable  to  the  Pearce  technique,  since  tins  error 
would  cancel  out.  This  is  one  advantage  of  using  amplitude  ratios  at  a  given  sta¬ 
tion  rather  than  the  amplitudes  themselves.)  In  order  to  speed  up  the  grid 
search,  we  have  also  revised  the  logic  so  that,  after  the  maximum  permitted 
number  of  inconsistencies  is  exceeded  for  a  given  grid  point  that  point  is  not 
included  in  any  further  searches  over  the  grid  This  "cascading"  of  acceptable 
solutions  from  one  ratio  test  to  another  made  the  grid  searching  several  times 
faster  Finally,  we  changed  Pearce's  treatment  of  the  polarity  data  to  consider 
the  polarities  of  the  phases  individually  rather  than  just  their  "ratio"  (i  e  , 
whether  the  two  phases  had  the  same  or  opposite  signs)  This  change  permits 
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more  P- wave  polarities  to  be  used,  since  usually  the  polarity  ...  'he  other  phase, 
whether  pP  or  sP,  cannot  be  observed,  and  therefore  no  test  can  be  performed 
on  the  "ratio"  of  the  polarities  of  the  two  phases. 

So  far  wc  have  described  three  techniques  for  finding  the  double-couple 
solution  (or,  more  precisely,  the  range  of  double-couple  solutions)  for  an  earth¬ 
quake  the  traditional  method  of  P- wave  (and  possibly  S’- wave)  polarities,  the 
use  of  amplitude  data  to  find  the  seismic  moment  tensor,  and  Pearce's  tech¬ 
nique  of  amplitude  ratios.  In  addition  to  their  use  for  finding  the  focal  mechan¬ 
isms  of  earthquakes,  each  of  these  techniques  may  also  be  used  as  a  tool  for 
discriminating  earthquakes  from  explosions  Since  our  principal  interest  in  this 
study  is  the  event  identification  problem,  we  shall  now  examine  the  use  of  each 
one  of  these  techniques  as  a  seismic  discriminant. 

The  traditional  technique  of  plotting  P- wave  polarities  on  a  sterconet  con¬ 
stitutes  a  simple  discriminant  :  if  an  event  has  at  least  one  dilatationa!  first 
motion,  then  it  is  an  earthquake.  This  straightforward  method  of  identifying 
earthquakes  was  in  fact  one  of  the  earliest,  discriminants,  having  been  proposed 
by  the  Geneva  "Conference  of  Experts"  in  19bB  An  appealing  aspect  of  this 
technique  is  that,  unlike  many  other  commonly  used  hut  empirically  derived 
discriminants  such  as  Ms . mb ,  its  physical  foundation  is  not  subject  to  contro¬ 
versy.  Unfortunately,  the  polarity  of  the  first  motion  for  weak  signals  is  usually 
ambiguous  and  sometimes  misleading  The  true  first  motion  can  be  obscured  by 
earth  noise,  causing  the  second  swing  of  the  seismograph  to  be  misread  as  being 
the  first  motion  For  this  reason,  the  first  motion  discriminant  is  applied  only  to 
those  waveforms  which  have  a  strong  signal-to-noisc  ratio  Consequently,  this 
technique  is  frequently  inadequate  for  identifying  small  cart : uptakes  Even  when 
the  signal-to-noisc  ratio  criterion  is  met,  the  discriminant  may  fail  since  the  /’ 
wave  polarities  of  certain  signals  from  explosions  nave  been  observed  to  be 
negative  (Knescu  et  al  ,  1973).  Ignoring  for  the  time  being  the  possibility  of 


misleading  first,  motions,  we  may  ask  whether  the  observation  of  only  compn-s- 
sional  polarities  for  an  event  means  that  the  event  is  an  explosion  In  most 
cases  the  answer  is  clearly  no,  since  a  pair  of  orthogonal  ares  representing 
nodal  planes  can  usually  be  drawn  on  the  stereonet  which  separate  the  stations 
into  a  pair  of  opposite  quadrants  in  which  compressional  polarities  are  observed 
and  a  pair  of  opposite  quadrants  in  which  no  polarities  are  reported  Only  if 
there  is  extensive  coverage  in  both  azimuth  and  distance  of  those  stations  which 
report  polarities  can  no  such  pair  of  nodal  planes  be  found.  If  it  does  turn  out 
through  a  complete  grid  search  of  all  possible  double-couple  mechanisms  that 
this  is  in  fact  the  case,  only  then  may  the  event  be  identified  as  an  explosion, 
otherwise,  an  event  with  only  compressional  first  motions  may  be  either  an 
explosion  or  an  earthquake. 

As  we  have  explained,  Pearce's  technique  uses  a  grid  search  over  all  possi¬ 
ble  double  couples  to  find  those  mechanisms  which  are  consistent  with  the 
observed  bounds  on  the  amplitude  ratios  and  polarities.  If  no  double  couple  is 
found  which  is  consistent  with  all  the  data,  then  we  conclude  that  either  some 
mistake  has  been  made  or  that  the  event  is  not  an  earthquake,  so  this  technique 
too  can  be  used  as  a  discriminant.  We  shall  discuss  in  a  subsequent  section  of 
this  report  possible  sources  of  error  which  could  cause  no  double-couple  solu¬ 
tion  to  be  found  for  an  earthquake.  For  the  moment  we  set  aside  the  question  of 
whether  an  earthquake  might  be  misidentified  by  Pearce's  technique,  and  ask 
whether  an  explosion  will  always  be  recognized  by  this  method.  In  general,  the 
answer  to  this  question  is  no,  since  in  analogy  to  the  previously  discussed  case  of 
a  stereonet  solution  which  involves  only  positive  first  motions,  we  see  that  usu¬ 
ally  there  exists  at  least  a  narrow  range  of  double-couple  solutions  which  arc 
consistent  with  all  the  data  unless  the  network  of  stations  provides  broad  azimu¬ 
thal  and  distance  coverage  Also,  in  analogy  to  the  case  of  small  events  for 
which  the  first  motion  is  usually  ambiguous,  for  a  small  explosion  pP  and  sP  will 
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be  obseured  by  earth  noise  The  "coda"  amplitude  which  i.  d  ,;s  an  upper 
bound  for  the  amplitude  of  the  undetected  surface  reflections  is  then  simply  the 
noise  level,  so  the  ratios  j>lJ/  I1  and  si3/  lJ  are  bounded  by  zero  ..rid  by  t  he 
reciprocal  of  the  signal-to-noiso  ratio  This  upper  bound  may  he  tro  large  to 
constrain  to  zero  the  range  of  consistent  solutions  In  this  study  we  shall  exam¬ 
ine  the  question  of  how  broad  the  coverage  of  stations  must  be  for  ! ’farce's 
technique  to  be  useful  as  a  discriminant,  and  we  shall  briefly  adores-  '  he  ques¬ 
tion  of  how  small  an  explosion  can  be  identified  by  this  method 

The  first-motion  criterion  has  been  used  extensively  as  a  seismic  discrim¬ 
inant  for  many  years  Inverting  seismic  measurements  to  solve  for  I  '~e  moment 
tensor  has  been  used  more  for  investigating  features  of  the  sc— -r  ■r.r'chamsn-.. 
such  as  the  yields,  of  known  explosions  than  for  event  identifF  1  ut  some 

studies  have  employed  this  technique  as  a  discriminant  (e.g  ,  S' r  ,'ev  -■ 

son,  1977)  Although  Pearce  has  applied  his  method  of  amplit1  !«•  *  bos  to  the 
problem  of  determining  the  focal  mechanisms  of  suites  of  cart'  ,  (Pc arcs 

et  al. ,  1980),  he  has  used  it  as  a  discriminant  for  only  a  few  spr ..  "  -  -  ,is  such 

as  the  20  March  1976  "Semipalatinsk  earthquake"  (Pooley  et  r.L. ,  '  9:T1)  The 
purpose  of  this  study  is  to  investigate  more  fully  the  utility,  from  both  a  theoret¬ 
ical  and  a  practical  standpoint,  of  this  technique  both  for  studying  earthquake 
focal  mechanisms  and  for  identifying  explosions.  In  particular,  .  e  shall  investi¬ 
gate  whether  the  necessary  measurements  for  implementing  Pearce's  algorithm 
may  be  made  in  a  manner  which  is  compatible  with  procedures  which  arc 
currently  used  for  routine  seismic  analysis. 

The  work  which  is  summarized  in  this  report  is  divided  into  throe  tasks 
First,  we  shall  present  the  results  of  an  investigation  of  certain  theoretical 
aspects  of  Pearce's  method,  such  as  the  azimuthal  coverage  which  is  r.cec.-sary 
to  constrain  tightly  the  range  of  solutions.  This  investigation  was  performed 
using  Global  Digital  Seismic  Network  [GDSN]  data  at  the  facilities  of  the  Center 
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for  Seismic  Studies  [CSS]  in  Arlington,  Virginia  The  use  of  the:  <  data  ai.d  f.a  .li 
ties  permitted  this  task  to  he  performed  in  an  unclassified  mode  Wxi  v..  pn. 
vide  documentation  of  a  package  of  computer  software  which  was  wr  - 

facilitate  making  the  measurements  on  seismograms  which  are  r<  •<,•<;;?<  d  h;, 
Pearce  s  algorithm  This  software  package  constitutes  a  new  subsystem  of  tin 
Seismic  Research  Information  System  [SHIS]  which  was,  at  the  time  of  tin--' 
study,  operating  in  a  classified  mode  at  the  Alexandria  Labs  of  Tt  ltd; 
Geotech  Finally,  we  present  results  of  a  demonstration  of  this  new  subsystem 
The  demonstration  was  performed  using  a  classified  data  base,  namely  the  Area 
of  Interest  [AI]  data  set  which  had  been  used  previously  in  the  YSC  Discrimina¬ 
tion  Experiment  and  which  was  used  concurrently  with  this  investigation  by 
another  study  involving  automated  discrimination  techniques  The  results  of 


this  classified  analysis  will  be  presented  in  a  separate  volume  of  this  report 
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The  seismologist  may  use  I  hi-  1  ‘ear  re  algonl  hm  ,  u  i-itiu-r  t>f  two  wavs  i' :  ‘.he 
first  wav,  a  seismologist  mu  v  simply  wish  to  tost  "the  earthquake  .it  depth 
hy  pot  ho  sis  of  .i  -aispeet  ed  event  This  tost  out  .ills  t  tie  est  uu.it  ion  of  I  ’  am:  i '  i  .vm 
love's  with  oa-eept  able  tolerances  The  amplitude  ratio  bounds  of  pi’  to  ami 
sI’lo-P  are  (lie  result  of  minimum  P  amplitude  bounds  arid  maximum  P  coda 
hounds  Simply  expressed,  the  pP-to-P  amplitude  ratio  bounds  are 


y/  coda  f  f  minimum  ^  -via/  ^  minimum) 

The  first  motion  polarity  of  the  P  wave  may  or  not  be  used  for  one  or  more  of  the 
stations  If  no  double-couples  can  be  found  consistent  with  the  pP/P,  -d1'!’ 
bounds  then  the  event  may  be  indicated  as  an  explosion  However,  if  double¬ 
couples  can  be  found  for  a  given  depth,  then  no  definite  conclusion  is  reached 
and  the  event  may  he  either  explosion  or  earthquake  it  may  be  necessary  to 
test  several  depth  intervals  It  is  not  possible  to  rule  out  the  nary  shallow  earth¬ 
quake. 


As  the  second  strategy,  the  seismologist  may  wish  to  find  possible  double- 
couple  solutions  at  a  given  depth  to  test  the  working  hypothesis  that  certain 
secondary  arrivals  arc  depth  phases.  There  may  be  several  alternative  combina¬ 
tions  of  secondary  arrivals  across  a  network  of  stations  resulting  from  coda 
fluctuations,  or  multiple  earthquake  sources  In  the  case  of  explosion  shot 
arrays,  there  may  exist  well  defined  secondary  arrivals  across  the  entire  array 
In  any  of  these  circumstances  the  seismologist  makes  the  necessary  measure¬ 
ments,  notes  the  inferred  depth  and  estimates  tolerances  on  the  amplitude 
measurement-:  11. eh  alternative  hypothesis  may  be  submitted  to  the  Pearce 
algorithm  The  processor  may  return  a  few  consistent  solutions,  numerous 
acceptable  solutions,  or  no  acceptable  solutions  The  results  may  eliminate 
some  of  the  competing  hypotheses 


In  the  following  examples,  various  strategies  are  illustrated  with  examples 
to  both  verify  that  the  Pearce  processor  yields  consistent  results  with  other 
focal  mechanism  methods  and  to  explore  some  of  the  assumptions  inherent  in 
the  algorithm. 

Figure  3A  shows  the  GDSN  P  wave  signals  of  an  Fast  Kazahk  explosion  of 
10/12/80,  event  1-19  in  Table  I.  Kaao  and  Zobo  are  complicated  wave  trams  due 
to  interference  of  multiple  arrivals  and  are  not  used  in  the  analysis.  P  wave 
amplitude  bounds  are  indicated  by  the  double  horizontal  bar  P  coda  levels  are 
also  indicated  by  the  single  horizontal  bars  These  amplitude  bounds  were  sub¬ 
mitted  to  the  Pearce  algorithm  as  P.  pP  and  sP  bounds  (  except  for  Grfo  and 
Gumo)  The  coda  amplitudes  were  multiplied  by  1.33  to  give  the  pP  and  sP 
bounds.  No  first  motion  polarity  information  was  used.  Xo  double-couple  solu¬ 
tions  were  found  to  be  consistent  with  the  20  |pP/P  and  sP/P,  amplitude  ratio 
bounds.  The  P  wave  takeofT  angles  for  each  station  are  shown  on  an  equal  area 
lower  hemisphere  projection  in  Figure  3B  (  P  velocity  of  6  km /sec  at  the 
source).  A  histogram  of  the  percentage  number  of  solutions  allowed  if  incon¬ 
sistencies  are  permitted  is  shown  in  Figure  3B.  For  example,  if  four  observa¬ 
tions  are  deleted,  then  96  solutions  out  of  '.1664  (  10  degree  search  grid  )  may¬ 
be  found.  This  presentation  of  the  results  does  not  show  which  observations  are 
critical.  However,  if  the  three  stations  Kono,  Anto,  and  Bcao  are  deleted  from 
the  data  set  then  an  180  degree  azimuthal  gap  opens  in  the  coverage  (Figure  3B) 
and  a  great  number  of  solutions  become  possible  The  P  coda  levels  are  valid  for 
a  depth  range  of  3  to  33  km.  No  acceptable  double-couple  solutions  were  found 
for  these  coda  levels  with  P  wave  velocities  of  b,  6  or  ?  km /sec  at  the  source  It 
is  quite  clear  that  for  this  event,  signal  to  ambient  noise  ratios  several  times 
worse  than  is  apparent  in  Figure  3A  would  still  discriminate  as  either  a  very 
shallow  (  less  than  3  km  deep)  earthquake  or  an  explosion  As  can  be  seen  from 


the  focal  sphere  projection  (Figure  3B),  that  stations  Anmo,  Lon  and  Jas  do  not 
contribute  independent  information  Similarly,  addition  of  Grfo  and  Gumo  only 
provides  redundant  information  to  Kono  and  Majo  respectively 

The  New  Brunswick  earthquake  serves  as  a  good  example  of  a  well  recorded 
continental  earthquake  with  published  solutions.  Figure  4 A  (  adapted  from  Ghoy 
et  al;  1983)  shows  the  short  period  GDSN  seismograms  and  broadband  decon¬ 
volved  records  used  by  Choy  et  al  (1983)  for  moment  tensor  inversion  of  the  New 
Brunswick  earthquake  Amplitudes  of  the  P  wave  and  depth  phases  pP  and  sP 
are  indicated  on  the  short  period  records  (SPZ)  as  they  were  used  in  the  Pearce 
focal  plane  analysis  program  Minimum  amplitudes  are  measured  for  pP  on 
ZOBO,  BOCO,  and  COL.  Minimum  amplitudes  are  measured  for  sP  on  TOL,  GRFO, 
and  KONO.  Maximum  amplitudes  are  given  for  pP  or  sP  on  each  SPZ  record  The 
polarities  were  considered  reliable  on  ZOBO.  and  BOCO  The  P  amplitude  bounds 
were  given  as  +  /-  50%  of  the  P  amplitude  measurements  and  submitted  as  input 
to  the  Pearce  algorithm.  Only  pP/P  and  sP/P  bounds  were  considered,  since  it 
was  not  clear  that  both  pP  and  sP  were  indentifiable  on  the  same  record  for  any 
of  the  stations  shown  in  Figure  4A.  Figure  4B  shows  the  results  of  a  5  degree 
search  in  slip,  dip,  and  strike  angles,  given  the  constraints  derived  from  the  SPZ 
records  of  Figure  4A.  64  focal  planes  (out  of  93312  possible)  were  found  con¬ 
sistent  with  any  9  out  of  the  12  observations.  When  pP/P  bounds  were  relaxed 
by  use  of  a  pP  attenuation  factor  of  0.75  slightly  expanded  set  of  solutions  was 
found  to  be  consistent  with  all  of  the  observations.  The  largest  number  of  solu¬ 
tions  were  found  for  a  P  velocity  at  the  source  of  5  km/sec  The  range  of  accept¬ 
able  solutions  can  be  described  as  having  a  steeply  dipping  tension  axis  and  a 
shallow  dipping  compression  axis.  The  azimuth  of  the  compression  axis  is  poorly 
constrained  by  the  data.  Figure  4C  shows  five  of  Choy  et  al's  solutions  for  this 
event,  and  the  Dziewonski  and  Woodhouse  (1963),  DW,  best  fitting  double-couple 
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axis.  Tension  (T)  and  Compression  (P)  axis  are  shown  on  a 
equal  area  lower  hemisphere  focal  sphere.  P  wave  take-off 
analyses  are  plotted  for  each  station  with  the  station  nam 
P  wave  nodal  planes  are  shown  for  each  acceptable  solution 
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These  published  solutions  <»ro  derived  from  several  cilh  'Methods  The  -solu¬ 
tions  suggest  that  the  bc*-t  interpretation  is  a  steeply  dipp  rig  ten,  -.r.  ax  is.  and  a 
shallow  dipping  I  last- West  compression  axis  Ail  but  the  ; )  solution  ire  con¬ 
tained  within  the  solution  space  shown  m  figure  Tne  !>W  .-oi.;;  t.;.  was  ba>r  d 
on  long  period  waveforms  and  was  admittedly  a  less  re  hub!  <  solution  The  C'f.oy 
et  ai  solutions  were  all  based  on  broad b  ind  body  waive*  data  To  sun.tv.a*  i/.e.  t  tie 
amplitudes  of  the  I*,  pi’  and  sP  phase-,:  ..pparent  on  six  Hi’/  records  at  Lelc-sc- 
isnne  distances  yielded  a  range  of  focal  mechanisms  very  dcse  to  that  of  much 
more  sophisticated  waveform  modeling  methods-  The  nerd  for  proper  attenua¬ 
tion  models  for  the  pP  and  sP  phases  :s  apparent  when  the  aim  m  to  find  the  best 
fitting  double-couple  solution  Alternatively,  the  error  bars  on  the  measured  P, 
pP,  and  sP  amplitudes  may  be  increased  The  experience  gained  with  this  exam¬ 
ple  was  that  increases  in  the  error  bounds  simply  increased  the  range  of  the 
acceptable  solution  sc!  for  a  g.ven  number  of  inconsistencies  tolerated  in  the 
observation  set  The  correct  attenuation  parameter  al'-owed  reduction  m  the 
error  bounds,  without  an  expansion  of  the  range  of  acceptable  solutions  The 
centroid  of  the  solution  space  was  relatively  stable  Tor  a  fixed  set  of  attenuation 
parameters,  while  the  volume  of  the  solution  spare  varied  with  the  size  of  the 
amplitude  error  bounds  that  were  used 

Two  events  of  special  interest  wih  recently  published  short  period  records 
are  the  events  nr  Vay  I  ,969  and  Varcn  .iO,  "976  near  .he  bast  ha/ahk  test  site 
r.smg  the  published  short  period  records  of  Pearce  et  ai  r  9dt'j  ,  and  Poo  ley  et  ,d 
we  wished  to  verify  that  the  previously  published  foeai  mechanisms  and 
depth-'  *•  (>>•(•  it.!  -.*.  ;,h  S  > , ,  •  ,iv..  'able  d  d  i  'Is  Tg  the  f  "ip!  o  -  l-'.ttinr'  of 

the  Pearce  dgorithm 

|'i-ic  pn  viousiv  pi.bi,  -nil  u  .mr,,-.  bi-.inh  I  '■*  '  h  A.  iS  a-\  a  no  ( IP  A  i  or  the  V  ay  ", 

;  0G9  event  are  shown  m  1'igun-  \  along  with  their  amplitudes  and  polarities 
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at  YKA 
t  =6.0° 
A=73.3° 


(b)  observed 
at  WRA 
i  =1276° 
A =81.6° 


(c)  observed 
at  GBA 
180.8° 
A=30.3° 
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FIGURE  1.  THREE  TELE SEISMIC  SHORT  PERIOD  ARRAY  OBSERVATIONS  OF  THE  1  MAY  1969 

EAST  KAZAKHSTAN  EA°THQUAKE,  ORIGIN  TIME  (K.00.08.7,  !,0CATI0N  43.98°N, 

77 . 86°E  AND  BODY  WAVE  MAGNITUDE  iil  =  U  .<)  (NEIS  PARAMETERS) .  T  AND  A* 
CORRESPOND  TO  AZIMUTH  OF  THE  STATION  FROM  THE  EARTHQUAKE  AND  EPICENTRAL 
DISTANCE  RESPECTIVELY.  PRESUMED  PHASE  IDENTIFICATIONS  AND  THEIR  POSSIBLE 
AMPLITUDES  ClN  ARBITRARY  UNITS)  ARE  SHOWN 


Figure  5a  Seismograms  reproduced  from  Pearce  (1979)  for  the  May  1,  1969 
event . 
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Figures  5B  and  5C  show  the  range  of  solutions  acceptable  with  these  amplitudes 
using  a  simple  one  layer  crustal  model.  The  north-west  striking  focal  planes 
may  be  compared  with  the  Chingiz  fault  which  strikes  N-40°-W  with  right 
lateral  offset.  The  fault  is  visible  on  LANDSAT  photos  reproduced  in  Rodean 
(1979).  Acceptable  double-couple  solutions  to  the  March  10.  1976  event  are 
shown  in  Figure  5D.  The  amplitudes  were  read  directly  from  recently  published 
seismograms  in  Pooley  et.  al.  (1983).  Pooley  et.  al.  arrive  at  similar  double¬ 
couple  solutions  for  the  March  20,  1976  event.  They  conclude,  however,  that  sP 
amplitudes  have  been  affected  by  lateral  heterogeneity  near  the  source  and  are 
less  reliably  fit  than  the  P  and  pP  amplitudes.  They  used  a  layered  crust  to 
predict  pP  and  sP  reflections  as  well  as  crustal  reverberation  in  their  simulation. 
Using  only  a  simple  one  layer  crust  without  reverberation,  the  synthetics  of  Fig¬ 
ure  5E.  reproduce  the  pP/P.  sP/P  and  pP/sP  amplitudes  ratios  observed  in  the 
actual  data  with  some  success.  In  genral  sP  amplitudes  are  at  the  upper  limit  of 
the  tolerances  submitted  to  the  Pearce  algorithm.  One  subset  of  the  bimodel 
solution  space  (Figure  5D)  supports  near  strike-slip  motion  upon  a  fault  striking 
NW-SE.  Pooley  et.  al. ’s  preferred  solutions  are  similar  to  those  in  Figure  5D 
Given  the  differences  between  a  simple  one  layer  crustal  model  and  a  layered 
reverberant  crustal  model  ,  this  exercise  shows  how  robust  the  Pearce  method 
may  be  to  overall  crustal  structure.  Discrepancies  are  most  likely  to  occur  with 
sP  amplitudes.  A  north-south  compressive  stress  system  would  be  consistent 
with  the  inferred  focal  mechanisms  of  both  the  May  1,  1969  and  March  20,  1976 
events. 

Two  earthquakes  near  30N,  95E  (Eastern  Tibet)  were  selected  with  short 
period  GDSN  recordings  (  events  1914  and  1716  in  Table  11).  Earthquakes  in  this 
area  were  previously  found  to  possess  explosion-like  M,  -m*  statistics  (Der, 
1973;  Landers,  1972).  Blandford  and  Gurski  (1975)  examined  available  geologic 
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•  Up  angla  in  fault  plana  PSI  (dag)  Plot  1914  • 

1914  95s29s82  11  27  22  8  29  669  95  469  33  8  89  8  99  8  •  l»eonaialant  oba 

19  39  69  99  129  159  189  out  of  8  oba 
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maps  and  1AMXSAT  image.-?,  and  hypothesised  that  these  anemalous  c.  c-rits  could 
be  small  ■'lb0  thrust  events  on  a  northwest-southca.i  striking  fault  system  with 
shallow  dip  The  available  seismic  data  has  been  parse  and  ir...  ufTicicnt  tr. 
determine  fault  planes  for  these  events  Seismogram-  of  the  l v. r<  events  art 
shown  in  figures  6A.M. C  Secondary  phases  in  the  coda  of  the  t wo  events  were 
assumed  to  be  pp,  sF5  arrivals  for  depths  of  \b  to  20  km  Results  of  the  Pearce 
algorithm  arc  presented  in  figures  61)  through  6i.  Tr-*  pf  *  P  and  sf.- 1’  amplitude 
ratio  bounds  were  found  to  be  consistent  with  thrust  sources  The  tension  axis 
was  well  constrained  to  the  near  vertical,  wniie  the  a,  :muth  of  the  near  horizon¬ 
tal  compression  axis  is  unconstrained  by  the  data  These  results  are  consist?  nt 
with  the  hypothesis  that  '".-V'  thrust  events  are  responsible  for  the  lew  /./4  meas¬ 
urements  of  the  TON  b>  9~>!f  events  Add  tior.ai  modeling  would  be  required  So 
test  if  the  compression;*!  axis  is  in  a  hen/ontai  nart he  as* -southwest  oriental  ion 
as  may  bo  inferred  from  lANPSAT  lir.em  mts 

Events  20  and  22  (label  II)  are  in  the  Take  Baikal  region  Event  20  has 
emergent  P  waves  and  is  possibly  a  double  event  (figui  os  ?A.H).  The  arrivals  .0 
seconds  following  the  P  wave  are  interpreted  as  pf  phases  The  best  fit  to  the 
data  is  a  horizontally  oriented  NE-.SW  tension  axis  with  n<  ariv  vertical  compres¬ 
sive  axis  (figure  7C).  The  event  is  200  km  west  of  lake  Baikal  and  apparent. h  *-ot. 
aligned  with  tne  rift,  zone  stress  system  It  was  not  passible  *<  find  double- 
couple  solutions  that,  fit  the  initial  two  f  phases  as  f  -pf  In  this  ease  the  focal 
mechanism  algorithm  was  used  to  decide  between  two  alternative  interpreta¬ 
tions  of  the  waveforms  The  event  appears  to  be  a  double  earthquake  with  «i mi- 
la*-  frx  •:  r.eeh.mi  -iii.'.  Th*-  sf  phase  i<  not  premia?  nt  i-:,  ‘  h**  •  >  a  .  -  and  rr'qt.ir?  d 
large  aniplit  tide  c  ops!  r:\int  s  The  simple  reflection  eoefTieir-nt  mode!  used  w..-- 
madequat?  to  mud*  I  these  sf  amplitudes  observed 

Event  22  (label  II)  is  located  within  the  Baikal  rift  zone  and  decidedly  shal- 
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lower  than  event  20  (Figure  8A).  A  solution  set  of  double-couples  with  horizontal 
NW-SE  tension  axis  and  near  vertical  compression  axis  was  found  consistent  with 
pP/P  and  sP/P  bounds  (Figure  8B).  Polarity  data  from  the  NW  azimuths  (Figure 
GC)  suggests  that  adjustments  may  be  necessary  if  waveform  data  is  added  in 
the  NW  quadrant  Although  the  focal  mechanisms  derived  from  pP/P  and  sP/P 
constraints  agree  with  the  expected  tectonic  stress  system,  it  may  not  satisfy 
additional  data  from  a  missing  quadrant 

Event  8  (Table  II)  occurred  on  the  arctic  coast  of  central  Siberia  in  an  area 
of  low  but  detectable  seismicity  (bang  and  Sun;  '.966)  Seismograms  and  accept¬ 
able  focal  mechanisms  based  on  pP  and  sP  bounds  are  shown  in  Figures  9A  and 
9B.  Focal  mechanisms  for  event  8  based  on  tight  pP  and  sP  bounds  are  con¬ 
sistent  with  nearly  all  the  polarity  data  of  Figure  9C.  Four  sP/P  amplitude 
bounds  must  be  violated  to  arrive  at  the  solutions  of  Figure  913 .  Scattered 
compressional  first  motions  suggest  that  the  east-west  striking  focal  planes  of 
Figure  9D,  might  be  rotated  to  a  more  nearly  ENE-WSW  strike  If  the  sP  bounds 
are  relaxed,  more  solutions  become  acceptable,  and  it  is  not  necessary  to  dis¬ 
card  data.  The  enlarged  solution  space  is  the  trade-ofT  for  Lhe  relaxed  sP 
bounds.  These  experiments,  show  that  the  Pearce  algorithm  may  be  used  in  an 
interactive  manner,  to  explore  the  robustness  of  the  results  In  this  case,  it  is 
clear  that  the  compressive  axis  is  steeply  dipping  to  the  northwest,  while  the 
tension  axis  is  not  as  well  constrained  If  the  ISC  polarity  data  were  added  to  the 
polarities  and  amplitudes  of  the  seismograms  of  Figure  bA,  the  preferred  ten¬ 
sion  axis  would  be  horizontal  NE-SW.  Although  the  location  of  this  event  is  far  to 
the  southeast  of  the  active  extension  of  the  Lomonosov  ridge,  the  stress  orienta¬ 
tion  appears  to  be  consistent  with  NE-SW  cxtonsional  tectonics  Also,  the 
results  indicate  that  the  sP  amplitudes  predicted  by  the  simple  crustal  model 
are  too  large  by  about  a  factor  of  2  for  this  source  depth  and  locality 


Event  295  (Table  II)  occurred  near  the  Kunar  fault  in  eastern  Afghanistan 
(seismograms  in  Figure  10A)  Prevot  et  al  (1980)  have  studied  the  tectonics  of 
this  area  and  found  focal  mechanisms  indicative  of  a  NW-SE  compressive  stress 
system.  A  NEIS  moment  tensor  solution  derived  from  long  period  waves  is  also 
available  for  this  event.  Figure  10B  shows  the  NEIS  double-couple  solution 
together  with  the  available  short  period  and  long  period  first  motions  reported 
to  NEIS.  The  NEIS  double-couple  solution  is  largely  in  agreement  with  the  first 
motions.  Crude  synthetics  for  the  NEIS  double-couple  solution  predict  large  pP 
and  sP  amplitudes  at  Nwao,  Bcao,  Grfo,  Kono,  and  Tol  (  Figure  10C)  Unless 
depth  phases  are  greatly  attenuated  in  this  source  region,  the  NEIS  double¬ 
couple  solution  is  inconsistent  with  the  P  coda  levels  observed  at  Nwao,  Bcao. 
Grfo,  Kono,  and  Tol.  Figures  10D  and  10E  are  the  acceptable  solutions  derived 
from  the  pP/P  constraints  on  the  traces  shown  in  Figure  10A  One  subset  of  the 
Pearce  solution  space  has  a  nearly  vertical  tension  axis  and  SE-NW  compression 
axis,  similar  to  the  NEIS  double-couple  solution  The  orientations  of  the 
compression  axes  are  significantly  different  enough  that  the  Pearce  solutions  do 
not  predict  large  pP  phases  The  Pearce  solutions  still  predict  sizable  sP  phases 
using  the  simplistic  crustal  structure.  The  discrepancy  is  at  least  a  factor  of  2. 
Allowance  for  acceptable  variations  of  the  S-to-P  velocity  ratio  at  the  source 
may  account  for  some  of  amplitude  discrepancy  However,  an  S  wave  velocity 
gradient  between  source  and  surface  will  reduce  the  SV  component  at  the  sur¬ 
face  reflection  point  dramatically,  and  hence  reduce  the  sP  amplitude. 

Event  516  (Table  II)  is  located  in  the  Tadzhik  region  and  has  an  an  indepen¬ 
dently  determined  double-couple  focal  mechanism  from  NEIS  The  GDSN  short 
period  seismograms  are  shown  in  Figure  1 1A  The  NEIS  reported  polarities  and 
double-couple  solution  are  plotted  in  Figure  1  IB  The  NEIS  double-couple  solu¬ 
tion  is  in  disagreement  with  most  of  the  reported  polarities  at  southerly 
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azimuths.  The  first  motions  are  largely  internally  consistent  and  indicate  that 
the  discrepancy  is  not  do  to  incorrectly  chosen  first  motions  Using  pi’  and  sF 
bounds  from  the  seismograms  of  Figure  11A.  the  acceptable  focal  mechanisms 
(Figure  1 1 C)  are  few  and  indicate  a  near  vertical  focal  plane  that  is  in  much 
better  agreement  with  the  polarity  data  (Figure  1  IB).  Also,  the  seismograms  of 
Figure  11A  indicate  that  Zobo,  Ctao,  Gumo,  and  Bcao  are  near  nodal  Again,  the 
NEIS  double-couple  solutions  derived  from  long  period  data  disagree  with  the 
short  period  data. 
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TKSTS  USING  SYNTIIhTIC  DATA 


In  order  to  bettor  understand  the  power  of  pl’/P  and  sI’/P  amplitude* 
bounds  as  constraints  on  the  allowable  range  of  focal  mechanisms,  several 
experiments  were  conducted  with  synthetic  data  These  experiments  were 
intended  to  explore  the  requirements  of  a  network  of  stations  to  eliminate  .ill 
possible  focal  mechanisms  in  the  case  of  an  explosion  In  practice,  maximum 
allowable  pi1,  and  si1  amplitudes  would  be  read  directly  form  the  coda  of  the 
explosion  I1  wave  and  would  be  limited  by  the  coda  generation  mechanisms  and 
the  ambient  noise  level 

Figures  12A.B  and  C  show  the  fraction  of  allowable  double-couple  solutions 
as  functions  of  relative  depth  phase  amplitude  and  take-off  angle,  for  various 
network  configurations  Figure  1 2 A  shows  the  fraction  of  allowable  focal 
mechanisms  versus  takcofT  angle  for  a  single  station  with  the  constraints  that 
|pPj  <  |P[,  |sP|  <  jPI,  and  both  ,pl’j,  and  ;sP  arc  <  'P;  It  is  apparent  that  tin*  take¬ 
off  angle  has  little  effect  for  the  case  of  jpl’/Pj  <  while  IsP/P,  <  I  is  sensitive  to 
takc-ofT  angle  However,  since  the  pP  constraint  is  more  powerful  than  the  sP 
constraint,  the  consideration  of  both  gives  little  reduction  in  the  possible  volu¬ 
tion  space  over  the  pP  constraint  alone  The  observation  that  depth  phase 
amplitudes  are  less  than  the  initial  phase  amplitude  at  a  single  station  elim¬ 
inates  roughly  50%  of  the  allowable  mechanisms  The  polarities  of  the  P  or  pP 
have  not  been  considered 

The  addition  of  stations  at  90  degree  azimuthal  intervals  is  addressed  in 
Figure  I2B  The  bounds  on  ]pP/P|  and  |sP/P|  were  allowed  to  vary,  and  the  frac¬ 
tion  of  allowable  solutions  was  plotted  for  the  cases  of  one,  two,  three,  and  four 
stations  at  90  degree  azimuthal  spacing  The  takc-ofT  angle  was  held  constant  at 
20  degrees.  The  polarities  of  the  P  waves  were  not  considered  We  see  that  if 
only  two  stations  are  observed,  then  the  coda  must  exhibit  levels  less  than  0  •'  of 


the  P  wave  to  eliminate  all  possible  double-couples  If  three  or  four  stations  are 
available,  then  coda  levels  of  0.95  or  less  are  sufficient  to  eliminate  all  possible 
double-couples 

The  results  of  a  similar  numerical  experiment  are  summarized  in  Figure 
12C,  for  a  network  of  stations  at  60  degree  azimuthal  intervals  and  20  degree 
take-ofT  angles.  The  results  are  somewhat  similar  to  the  previous  experiment,  in 
that  four  or  more  stations  with  coda  levels  slightly  less  than  the  P  wave  are 
sufficient  to  eliminate  all  possible  double-couple  solutions.  The  differences 
between  the  two  results  indicate  that  the  network’s  azimuthal  coverage  may  be 
important  to  the  inherent  capability  of  the  network  to  discriminate  between 
earthquakes  and  explosions.  To  simulate  something  closer  to  the  azimuthal  cov¬ 
erage  that  may  be  available  in  practice,  additional  simulations  were  conducted 
with  a  random  distribution  of  azimuths  and  teleseismic  take-off  angles.  Further¬ 
more,  since  the  coda  levels  experienced  in  practice  are  not  the  same  at  all  sta¬ 
tions,  a  random  model  of  the  |pP/P|  levels  was  used. 

A  Monte  Carlo  simulation  of  a  random  network  with  random  pP/P  amplitude 
ratios  was  performed.  Five  or  more  realizations  of  20  stations  was  generated  for 
each  set  of  statistics.  The  station  azimuths  were  evenly  distributed  in  360 
degrees  of  azimuth  and  take-off  angles  were  evenly  distributed  between  10  and 
35  degrees.  The  |pP/P|  amplitude  ratio  bounds  were  given  by  a  normal  distribu¬ 
tion  with  means  of  1.0,  0.75  and  0.5,  and  standard  deviations  of  0.25  and  0.5.  The 
number  of  stations  required  to  eliminate  all  double-couples  was  counted  for 
each  realization  A  survey  of  the  simulations  indicates  that  explosions  will  be 
inconsistent  with  double-couples  when  three  out  of  four  azimuthal  quadrants  are 
represented  with  coda  amplitude  levels  below  the  P  amplitude  level  Since  270 
degree  azimuthal  coverage  will  require  four  or  more  stations,  the  simulations 
came  down  to  the  probability  that  four  stations  would  be  "drawn  out  of  the  hat’’ 
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Figure  12a  The  percentage  of  all  possible  double-couples  eliminated  by 
the  constraint  of  |pP/p|  <  1 ,  and  |sP/P|  <1  at  a  single 
station.  Constraints  are  considered  together  and 
separately.  The  teleseismic  takeoff  angle  is  varied  from  10 
to  25  degrees.  The  observation  that  |  pP  /  P |  <  1  eliminates 
nearly  507  of  all  double-couples  at  any  teleseismic  takeoff 
angle.  The  implications  that  sP  p  are  less  significant  and 
slightly  dependent  on  takeoff  angle.  The  observation  that 
both  pP  and  sP  are  smaller  than  P  implies  that  slightly 
fewer  solutions  are  acceptable  the  observation  the  pP  <P 
alone. 


*/.  of  focal  mechanisms  allowed 


90  DEGREE  AZIMUTHAL  INCREMENTS 


0.0  1.0  1.5 

pP/P,  sP/P  limits 


Figure  12b  The  percentage  of  all  possible  double-couples  eliminated  by 
one,  two,  three,  or  four  stations  at  90  degree  azimuthal 
increments  varies  as  a  function  of  the  absolute  limit  of  the 
pP  and  sP  amplitudes.  All  takeoff  angles  are  assumed  to  be 
20  degrees. 


60  DEGREE  AZIMUTHAL  INCREMENTS 


pP/P,  sP/P  limits 


Figure  12c  The  percentage  of  all  possible  rloubl  e- couples  eliminated  bv 
a  network  of  stations  at  60  degree  azimuthal  increments. 

Four  or  more  stations  are  required  with  |  pP/P  |  and  |  sP/P|^  1 
to  eliminate  all  possible  focal  mechanisms  from 
consideration. 


with  sufficiently  small  coda  amplitudes.  The  stochastic  variation  in  teleseismic 
take-off  angles  did  not  seem  to  be  a  significant  factor.  The  average  number  of 
stations  required  to  eliminate  all  double-couples  for  a  mean  |pP/P|  of  0.75  with 
0.25  standard  deviation  was  8  stations.  The  average  number  of  stations  required 
to  eliminate  all  double-couples  for  a  mean  |pP/P|  of  0.5  with  0.25  standard  devia¬ 
tion  was  6  stations.  For  jpP/Pj  =  0.75  with  a  standard  deviation  of  0.5  an  average 
of  11  stations  was  required  and  the  results  were  highly  variable.  Values  of  |pP/P| 
>  1.0  were  of  little  use  in  eliminating  double-couple  solutions. 

As  one  last  set  of  tests,  a  series  of  hypothetical  shot  arrays  were  investi¬ 
gated  with  a  nine  station  network  appropriate  to  the  east  Kazahk  test  site.  Four 
hypothetical  shot  ratios  were  considered;  the  first  shot  equal  in  magnitude  to 
the  second  shot,  ie  shot l/shot2= 1.0.  as  well  as  the  ratios  shotl/shot2=0.5, 
shotl/shot2=.66,  shotl/shot2=1.5.  and  shotl /shot2=2.0  The  extremal  bounds 
on  the  two  arrivals  were  assumed  to  be  +/-  25  percent.  The  only  shot  ratio  that 
allowed  solutions  was  shotl/shot2=l.  The  nine  station  test  array  together  with 
the  solution  set  for  (l  +  /-0.25)/(l  +  /-0.25)  is  shown  in  Figures  13A  and  B.  These 
solutions  would  be  easily  removed  if  polarity  information  was  available  on  the 
first  shot  since  the  double-couple  solutions  all  possess  vertical  nodal  planes,  and 
would  be  eliminated  by  any  three  stations  with  the  same  polarity  separated  by 
90  degrees  of  azimuth.  Results  for  the  other  shot  ratios  are  robust  with  respect 
to  deletion  of  one  or  two  stations  from  the  network.  The  Pearce  algorithm  would 
appear  to  have  some  potential  for  discrimination  of  shot  arrays  designed  to 
simulate  P+pP  if  the  two  shots  were  of  greatly  different  magnitude,  or  if  the  two 
shots  were  of  comparable  magnitude  and  polarities  were  available  for  two  or 
more  quadrants. 
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Figure  13a  Results  of  a  hypothetical  nine  station  network  with  two 
arrivals  of  nearly  equal  amplitude.  Both  P  and  the 
prospective  pP  are  given  25%  amplitude  tolerances.  4626 
possible  focal  planes  out  of  93312  are  found.  No  polarity 
information  is  used.  Solutions  are  either  of  dip-slip  with 
normal  thrust,  dip-slip  with  reverse  thrust,  vertical 
dip-slip,  or  horizontal  faults  with  horizontal  slip. 
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DATA  PROCESSING  AT  THE  CSS 


Procedures  and  software  were  developed  at  the  CSS  to  input  amplitude 
measurements  of  P  arrivals  and  prospective  depth  phases  into  the  Pearce  algo¬ 
rithm  and  to  analyze  the  results  in  graphical  form.  Software  documentation  in 
the  form  of  manual  pages  is  reproduced  in  the  appendix.  A  flow  chart  of  the  typ¬ 
ical  seismic  analysis  is  shown  in  Figure  14.  The  hypocentral  location  may  be 
preliminary,  or  poorly  constrained  (33  km).  The  process  of  depth  estimation 
and  depth  phase  identification  may  be  an  interactive  process  where  the  seismol¬ 
ogist  1)  selects  prospective  depth  phases,  2)  infers  a  depth,  3)  estimates  the 
amplitude  bounds  of  P,  pP,  and  sP,  4)  runs  the  Pearce  algorithm,  and  5) 
analyzes  the  output.  The  seismologist  may  decide  that  either  the  depth  phases 
Eire  incorrectly  selected,  that  the  amplitude  estimates  may  be  improved,  or  that 
the  results  are  satisfactory. 

For  the  purpose  of  discrimination,  amplitudes  are  read  from  the  network  of 
stations  at  selected  time  windows  following  the  P  onset.  Each  Lime  window 
implies  a  depth  range,  and  the  Pearce  processor  is  run  as  a  batch  job  for  each 
depth  interval  in  consideration.  The  output  was  analyzed  at  a  later  time 

Although  the  UNIX  operating  environment  at  the  CSS  was  satisfactory  for 
the  interactive  aspects  of  focal  plane  analysis,  such  as  amplitude  measure¬ 
ments,  or  graphics  display,  the  Pearce  processor  was  too  slow  to  run  in  real 
time.  The  processor  was  typically  used  in  a  "quasi-batch"  mode  as  an  overnight 
process.  Restructuring  of  the  Pearce  algorithm  may  reduce  computation 
slightly.  However,  the  major  computational  task  of  calculating  nearly  100,000 
radiation  patterns  for  each  observation  (  5  degree  search  increments  ),  is  cen¬ 
tral  to  the  algorithm  and  not  easily  avoided.  An  option  to  reduce  computation 
was  added  to  eliminate  consideration  of  solutions  that  have  already  been  found 
inconsistent  with  more  than  a  set  number  of  observations.  This  option  reduces 


execution  time  at  the  expense  of  complete  statistics  for  each  possible  solution 

For  the  purpose  of  this  preliminary  study,  the  individual  seismograms  were 
interpreted  with  an  interactive  graphics  display  program.  The  amplitude  meas¬ 
urements  were  made  manually  by  control  of  a  cross-hair  terminal-host  com¬ 
puter  interface  In  the  case  of  most  explosions  and  a  few  earthquakes  where 
prospective  depth  phases  are  not  obvious,  an  automated  processor  could  be 
developed  that  gives  the  bounds  on  pP/P  for  given  set  of  depth  intervals  at  each 
station.  The  output  would  then  reviewed  by  the  seismologist  before  submission 
to  the  Pearce  processor  Such  an  automated  procedure  would  speed-up  the  pro¬ 
cedure  as  a  discrimination  tool.  If  such  a  processor  is  developed,  it  should  allow 
for  uncertainties  in  the  P  arrival  time,  and  use  both  seismologist  provided  P 
picks  as  well  as  predicted  arrival  times  based  on  hypocentral  parameters. 

Two  waveform  data  bases  were  used  to  evaluate  the  Pearce  focal  plane 
analysis  algorithm  at  the  Center  for  seismic  studies.  The  first  data  base  con¬ 
sisted  of  explosions  recorded  at  SRO,  ASRO,  and  DWWSSN  stations  The  under¬ 
ground  explosions  were  located  in  East  Kazahk,  Novalya  Zmelya,  and  Tuamoto 
test  sites.  The  majority  were  located  in  the  East  Kazahk  test  sites  The  second 
data  base  consisted  o'  Eurasian  earthquakes  recorded  on  the  same  network  of 
stations  available  from  GDSN  day  tapes.  Pertinent  information  for  these  two 
sets  of  events  are  presented  in  Tables  1  and  II.  The  waveforms  for  each  event 
were  displayed  and  measurements  of  the  P  wave  and  P  wave  coda  were  made 
with  two  objectives.  The  first  objective  was  to  test  if  a  simple  measure  of  the  P 
coda  maximum  would  discriminate  explosions  form  earthquakes  when  normal¬ 
ized  to  the  P  wave  amplitude  and  used  as  a  maximum  bound  on  p P/P  amplitude 
ratios  across  the  network  of  stations.  P  coda  was  defined  for  the  purpose  of 
these  measurements  as  that  portion  of  the  seismogram  more  than  3  seconds 
after  the  P  wave  onset.  This  corresponds  to  minimum  pP-P  times  of  3  seconds 

-79- 

*  1 1  ■  1  -  mm - * 


if 


or  dcplhs  of  about  b  km  or  greater  Since  the  time  of  the  prospective  pi* 
arrivals  on  each  seismogram  were  not  necessarily  consistent  with  some  fixed 
depth,  this  simplistic  estimate  of  the  maximum  |pP/P,  and  jsl’/l’  ratios  is  con¬ 
sidered  conservative.  Because  more  stringent  bounds  on  ipP/P  and  sP/P 
ratios  are  derived  from  forcing  some  consistency  with  a  fixed  depth,  a  second 
set  of  experiments  were  performed  to  find  the  best  estimates  of  depth  and 
source  mechanism  that  satisfied  the  prospective  depth  phase  arrivals  of  the 
earthquakes  and  explosion  waveforms.  For  the  earthquakes  and  explosions,  a 
Lime  window  following  the  P  wave  was  selected  with  a  large  excursion  in  the  coda 


on  one  or  more  seismograms.  Each  waveform  amplitude  was  measured  within  5 
seconds  of  the  prospective  depth  phase  arrival  time.  The  focal  plane  solutions 
for  the  earthquakes  were  compared  with  available  published  source  mechanisms 
in  the  the  same  region  and  with  polarity  data  reported  to  ISC.  and  NEIS. 
Selected  earthquakes  are  discussed  in  the  section  on  SOME  DETAILED  EXAM¬ 


PLES  OF  THE  PEARCE  PROCESSOR 

Ideally,  each  event  should  be  recorded  at  five  or  more  stations  at  distances 
between  30  and  90  degrees  with  uniform  azimuthal  coverage  Given  such  a 
situation  it  would  be  possible  to  test  the  hypothesis  that  each  event  was  an 
earthquake.  Stations  clustered  at  specific  azimuths  from  the  epicentral  loca¬ 
tion  do  not  give  independent  information  and  may  give  redundant  information 
Because  of  the  non-uniform  teleseismic  coverage  of  the  Eurasian  continent  by 
the  SRO,  ASRO  network,  and  numerous  down-time  exhibited  by  these  stations, 
the  addition  of  DWWSSN  stations  to  the  digital  network  day  tape  in  19B2  strongly 
enhanced  the  data  coverage.  Regional  distances  were  not  used  because  depth 
phase  arrival  times  are  more  sensitive  to  the  regional  velocity  structure  and 
because  the  upper  mantle  P  wave  travel  time  triplications  complicate  the  selec¬ 
tion  of  depth  phases.  For  the  same  reason,  some  PkP  observations  were  not 
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used  when  PkP  triplications  produced  interfering  secondary  arrivals  Available 
CSS  software  was  used  to  display  the  P+Pcoda  waveforms  with  predicted  P,  pP, 
sP,  PcP.  and  PkP(df)(bc)(ab)  arrival  times.  The  software  was  not  capable  of 
predicting  pPcP,  sPcP,  pPkP,  or  sPkP  arrival  times.  Nor  were  the  travel  Lime 
tables  adequate  for  shallow  events  (  less  than  33  km).  Differential  sP-pP  times 
were  usually  incorrect  for  shallow  events  and  the  pP  times  systematically  late 
Because  hypocentral  coordinates  were  usually  preliminary  and  because  of  local 
station  travel  time  biases,  the  predicted  P  and  depth  phase  times  were  often 
incorrect.  The  available  software  did  not  permit  an  interactive  relocation. 

As  a  first  pass  experiment,  explosions  and  earthquakes  were  tested  with 
depth  phase  amplitude  bounds  of  the  form;  |pP|  <  C.  The  value  of  C  was  read 
from  the  maximum  P  coda  amplitudes  following  the  P  arrival.  The  explosion 
amplitudes  were  submitted  to  the  Pearce  processor  for  takeoff  angles  with 
source  velocities  of  5.0,  and  6.0  km/sec.  The  earthquake  amplitudes  were  sub¬ 
mitted  to  the  processor  with  velocities  appropriate  for  the  best  guess  depths. 
None  of  the  earthquakes  failed  to  yield  some  solution  set  of  double-couples.  Out 
of  the  41  explosions  in  Table  I,  12  were  consistent  with  double-couple  solutions  in 
the  crust.  Only  two  of  these  12  events  had  5  or  more  stations.  All  of  the  other 
10  events  were  only  recorded  on  four  or  less  stations.  All  of  these  remaining  10 
events  were  not  recorded  in  one  or  more  quadrants.  Six  of  these  10  events  were 
recorded  in  only  2  quadrants  (  stations  located  only  in  Europe  and  North  Amer¬ 
ica  ).  Of  the  two  events  with  5  or  more  stations  recorded;  event  72  was  located 
in  the  Tuamato  islands  with  recordings  in  only  2  quadrants,  and  event  115  was 
also  recorded  in  only  2  quadrants.  These  12  explosions  spanned  the  magnitude 
range  of  5.2  to  6.0  m$.  Bandpass  filtering  was  used  on  seismograms  from  10  of 
the  41  explosions.  With  one  or  two  exceptions,  bandpass  Altering  did  not  reduce 
signiffcantly  depth  phase  amplitude  constraints.  The  early  P  coda  did  not  pos- 


sess  significantly  different  frequency  bandwidth  than  the  F  wave,  and  only  where 
ambient  noise  was  high  did  filtering  significantly  affect  the  amplitude  bounds. 

Several  of  the  earthquakes  with  well  constrained  focal  mechanisms  arc  dis¬ 
cussed  in  a  previous  section.  The  minimum  amplitude  bounds  for  pP  and  sP 
were  considered  as  well  as  the  maximum  amplitude  bounds,  P  polarities,  and 
more  stringent  bounds  on  depth  phases  were  interactively  varied  to  attain  a 
best  estimate  of  the  solution  set  of  focal  mechanisms.  Those  of  special  interest, 
with  independent  estimates  of  the  focal  mechanism,  or  illustrated  some  aspects 
of  algorithm  are  treated  with  detail  in  a  previous  section. 


PEARCE  ANALYSIS  SUBSYSTEM  OF  THE  SKIS 


In  the  previous  section  of  this  report  we  described  a  system  for  performing 
Pearce's  analysis  using  the  VAX  1 1/7B0  computers  which  run  under  the  UNIX  operating 
system  at  the  CSS.  In  this  section  we  shall  describe  another  system  for  performing 
Pearce’s  analysis  which  has  important  operational  differences  from  the  system  at  the 
CSS.  This  system  also  uses  a  VAX  11/780  computer,  but  it  runs  at  a  different  facility, 
the  Seismic  Research  Center  [SRC]  in  Alexandria,  Virginia,  and  it  uses  a  different 
operating  system,  namely  VMS  [Virtual  Memory  System].  The  system  at  the  SRC  is  a 
subsystem  of  a  larger  seismic  analysis  system,  whereas  the  CSS  system  operates  in  a 
stand-alone  mode.  The  two  systems  also  differ  in  two  respects  with  regard  to  their 
databases.  At  the  time  of  this  study,  the  facilities  at  the  CSS  were  restricted  to  the  use 
of  only  unclassified  data,  but  the  SRC  system  was  run  in  a  classified  mode.  Also,  the 
CSS  system  uses  waveforms  copied  to  disk  from  the  GDSN  tapes,  whereas  the  SRC  sys¬ 
tem  employs  a  highly  structured  file  system  consisting  of  event  records  and  waveform 
files  which  are  accessed  using  the  Record  Management  System  [RMS]  "indexed  file"  I/O 
features  of  VMS. 

The  contractual  task  which  was  fulfilled  by  the  construction  of  this  second  system 
was  that  of  adding  the  capability  of  performing  Pearce's  analysis  to  the  Seismic 
Research  Information  System  [SRIS]  (formerly  known  as  the  Regional  Event  Location 
System  [RELS]),  a  large  software  package  for  seismic  analysis  which  runs  on  the  VAX  at 
the  SRC.  Documentation  of  SRIS  has  been  published  elsewhere  (Teledyne  Geotech, 
1982),  and  we  shall  herein  assume  familiarity  on  the  part  of  the  reader  with  the  design 
and  operation  of  this  system.  In  its  original  formulation.  SRIS  contained  one  subsys¬ 
tem  devoted  to  the  discrimination  task,  namely  the  Discrimination  and  Identification  of 
Seismic  Events  [D1SE]  subsystem  (von  Seggern,  1981),  commonly  known  as  the 
"interactive"  discrimination  subsystem.  The  performance  of  that  subsystem  has  been 
tested  by  Blandford  et  oi.  (1983)  using  measurements  which  were  made  on  continuous 
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waveforms  by  a  new  SHIS  subsystem,  the  Automatic  Seismic  Discrimination  System 
[ASDIS]  (Farrell,  1981),  In  the  final  formulation  of  SKIS,  this  new  ''automatic"  discrimi¬ 
nation  subsystem  was  extended  to  handle  the  segmented  waveforms  of  the  A!  data  set, 
and  a  third  discrimination  subsystem  was  added  so  that  Pearce's  technique  could  be 
implemented  on  the  system.  We  shall  refer  to  this  last  subsystem  simply  as  the 
"Pearce"  subsystem.  In  order  to  provide  results  which  could  be  compared  with  previ¬ 
ous  studies  and  with  the  concurrent  testing  of  the  automatic  discrimination  subsys¬ 
tem.  the  Pearce  subsystem  was  designed  to  operate  on  segmented  data  files  so  that  it 
too  could  be  applied  to  the  A1  data  set.  In  its  final  configuration,  then,  SKIS  contained 
three  discrimination  subsystems  :  an  interactive  discrimination  subsystem,  which 
operated  on  alphanumeric  data  files;  an  automatic  discrimination  subsystem,  which 
originally  operated  only  on  continuous  waveform  data  but  later  on  segmented 
waveform  files  as  well,  and  which  generated  the  input  for  the  interactive  discrimination 
subsystem;  and  the  Pearce  discrimination  subsystem,  which  operated  on  segmented 
waveform  files.  (As  we  shall  discuss,  provision  was  also  made  for  the  Pearce  subsystem 
to  access  continuous  waveform  data,  but  this  feature  of  the  design  was  not  imple¬ 
mented.)  In  this  section  we  shall  present  certain  documentation  pertaining  to  the 
design  and  operation  of  the  SRIS  Pearce  discrimination  subsystem. 

The  Pearce  analysis  at  the  SRC  consists  of  three  separate  stages  of  data  process¬ 
ing.  First,  the  interactive  graphics  capabilities  of  the  SRIS  are  used  to  access  the 
desired  waveforms  and  to  make  the  amplitude  and  polarity  measurements  which  are 
needed  by  the  Pearce  algorithm.  Next,  this  algorithm  itself  is  run  as  a  batch  job  on 
VMS.  As  we  have  described  in  the  previous  section,  the  grid  searching  which  is  per¬ 
formed  by  this  algorithm  can  be  a  time-consuming  process,  particularly  when  the 
focal-mechanism  parameters  are  searched  in  increments  of  5°.  Cascading  the  grid 
search  from  one  amplitude  test  to  another  in  the  manner  which  we  have  described 
speeds  up  this  process,  but  it  also  loses  information  about  which  stations  most  strongly 
restrict  the  range  of  allowed  solutions  (this  information  may  be  valuable  for  identifying 
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critical,  or  perhaps  erroneous,  measurements),  and  even  with  "cascading"  it  is  still 
slow  enough  that  running  it  on-line  would  cause  a  troublesome  delay  before  the  SHIS 
would  be  free  again  for  further  waveform  processing.  Making  the  Pearce  algorithm  a 
batch  job  gives  the  SR1S  operator  the  option  of  either  deferring  this  analysis  until 
sometime  after  he  has  finished  using  SRIS  (in  which  case  he  may  submit  the  job  from  a 
"dumb"  terminal  not  connected  to  SRIS)  or  of  submitting  the  job  interactively  from  the 
SRIS  keyboard  and  monitoring  its  output  on  the  raster  screen  while  it  runs.  We  antici¬ 
pate  that  the  operator  will  select  this  second  option  if  the  Pearce  analysis  is  the  only 
processing  which  remains  to  be  performed  during  a  given  SRIS  session;  in  this  case,  the 
operator  will  want  to  know  while  he  is  still  at  the  SRIS  analysis  station  whether  the 
Pearce  algorithm  was  unsuccessful  and  whether  he  should  repeat  any  of  the  waveform 
measurements.  Finally,  the  range  of  solutions  allowed  by  the  Pearce  algorithm,  which 
is  stored  on  the  system  as  a  binary  file  so  that  it  can  be  accessed  at  any  time  in  the 
future,  is  input  to  an  interactive  job  which  displays  these  solutions  graphically.  Two 
versions  of  this  graphics  job  exist  :  one  uses  the  Tektronix  storage  tube,  which  is  not 
connected  with  SRIS.  and  the  other  uses  the  Graphicus-80  [G-80]  refresh  screen  which 
is  a  part  of  the  SRIS  analysis  station.  The  Tektronix  version  would  be  used  for  ofT-line 
processing,  such  as  for  displaying  the  results  of  the  Pearce  analysis  which  had  been 
run  as  a  batch  job  or  for  re-displaying  the  results  of  analysis  which  had  been  performed 
sometime  in  the  past.  The  G-80  version  would  be  used  for  on-line  processing,  so  that 
the  SRIS  operator  could  examine  the  results  of  a  Pearce  algorithm  job  which  he  had 
just  finished  running. 

We  shall  now  describe  the  features  of  the  new  SRIS  subsystem  which  permits  the 
operator  to  make  the  waveform  measurements  which  are  to  be  used  by  the  Pearce 
algorithm.  This  new  subsystem  is  an  extension  of  the  Interactive  Analysis  Station  [IAS] 
subsystem,  which  is  described  in  the  SRIS  documentation  volumes  (Teledyne  Geotech, 
1982).  This  new  subsystem  is  invoked  by  requesting  the  SUBMENU  feature  from  among 
the  other  options  which  are  offered  to  the  operator  on  the  "main  menu”  of  the  IAS, 
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Option 

GETCONWF 

GETSEGWF 

AMPUM 

POLARITY 

SCRSAVE 

T-P1CK 

DEPTH 

TAKEOFF 

DISPLAY 

P-P1CK 

MAINMENU 

IDLE 


SUBMENU  Options  of  the  Pearce  Measurement  Subsystem 


Function  Performed 

Retrieve  continuous  waveform  data 

Retrieve  segmented  waveform  data 

Make  measurements  of  maximum  or  minimum 
bounds  on  the  amplitude  of  signals  (or  of  coda) 

Make  picks  of  signal  polarity  or  0) 

Save  results  measured  so  far  into  scratch  file 

Display  time  corresponding  to  cursor  position  on 
waveform 

Choose  a  trial  value  of  the  event  depth;  for  this  depth, 
mark  the  predicted  arrivals  of  pP,  sP,  and  PcP 

Compute  the  takeoff  angle  of  the  P  phase  of  the 
waveform  being  measured 

Display  the  retrieved  waveforms  and  mark  the 
predicted  P  arrivals 

Replace  the  predicted  P  marker  with  an  observed 
pick 

Return  to  the  main  menu  of  the  IAS  to  perform 
waveform  processing 

Submenu  element  which  cannot  be  picked  but  which 
blinks  as  a  prompt  whenever  the  system  is  waiting  for 
a  submenu  pick 


which  is  displayed  on  the  G-80  screen.  Requesting  the  SUBMENU  clears  the  main  menu 
from  the  screen  and  replaces  it  with  a  new  menu  of  options,  all  of  which  pertain  to  the 
Pearce  measurements.  These  options  are  listed  and  defined  in  Table  111.  Strictly 
speaking,  the  IDLE  element  is  not  a  submenu  option;  it  merely  blinks  when  the  system 
is  in  a  state  in  which  it  can  accept  and  process  a  submenu  pick.  We  shall  now  describe 
how  an  analyst  can  use  these  options  to  make  the  waveform  measurements  which  are 
required  by  the  Pearce  algorithm. 

The  operator  will  first  need  to  access  the  waveform  data  pertaining  to  the  event 
under  analysis.  He  may  retrieve  continuous  waveform  data  using  submenu  option 
GETCONWF  or  segmented  data  using  option  GETSEGWF.  The  GETCONWF  option,  how¬ 
ever,  was  never  implemented  in  this  design,  since  we  were  tasked  to  demonstrate  the 
system  using  the  AI  (segmented)  dataset  rather  than  continuous  waveform  data.  This 
menu  option  was  thus  included  only  as  a  stub,  so  that  it  could  be  added  to  the  system 
at  a  later  time  if  it  was  desired  to  do  so.  This  option,  if  it  were  to  be  implemented, 
would  use  much  of  the  same  software  for  accessing  the  continuous  waveform  files  which 
is  used  by  the  IAS  and  by  the  on-line  Detection  Processing  [DP]  SR1S  subsystem.  That 
this  existing  software  may  in  fact  be  used  successfully  for  retrieving  continuous 
waveform  data  has  been  demonstrated  by  Blandford  et  a l.  (1983).  We  shall  therefore 
desribe  only  the  option  GETSEGWF,  which  employs  a  newly  developed  software  package 
for  retrieving  segmented  waveform  files  such  as  those  of  the  AI  dataset. 

The  GETSEGWF  option  employs  two  types  of  data  files  which  are  described  in  the 
Data  Base  Specification  [DS]  volume  of  the  SRIS  documentation  (Teledyne  Geotech, 
1982)  :  Event  Arrival  Files  [EAF]  and  Segmented  Data  Files  [SDF].  The  operator  directs 
the  system  to  read  the  proper  EAF  record  by  replying  to  its  prompt,  which  appears  on 
the  raster  screen,  to  enter  the  year,  Julian  day  within  the  year,  and  event  number 
within  the  day  of  the  desired  event.  The  event  arrivals  for  each  day  constitute  a  single 
EAF,  so  the  first  two  items  entered  are  sufficient  to  cause  the  system  to  open  the 
proper  VMS/RMS  "indexed"  file.  (It  is  also  necessary  for  the  system  to  know  on  what 
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disk  and  within  what  directory  the  EAF  resides;  this  information  is  hard-coded  within 
subroutine  RDSEG.)  All  three  of  the  items  entered  are  concatenated  to  form  the  "key" 
which  enables  the  RMS  feature  of  VMS  to  find  the  correct  record  in  the  indexed  file 
The  system  reads  the  "type  0"  EAF  record  for  the  sought-after  event,  which  includes 
such  information  as  the  event  hypocenter  and  origin  time.  Information  about  the  indi¬ 
vidual  signal  arrivals  from  this  event  as  recorded  at  stations  throughout  the  network 
are  contained  in  the  "type  1"  records  which  follow  this  type  0  record  in  the  EAF.  The 
system  reads  through  these  records  (by  reading  the  first  one  with  a  keyed  read  and 
then  using  sequential  reads  after  that)  in  order  to  find  out  what  waveforms  are  avail¬ 
able.  Before  reading  the  EAF  type  1  records,  the  system  prompts  for  the  desired  phase 
(for  our  purposes.  P),  and  it  opens  the  SDF  for  the  same  Julian  day  as  that  of  the  EAF. 
Consequently,  when  it  finds  an  EAF  type  1  record  of  a  phase  of  the  designated  type 
(i.e.  .  short-period  P),  it  can  perform  a  keyed  read  of  the  corresponding  SDF  record  to 
retrieve  the  waveform  containing  that  arrival.  The  EAF  thus  functions  as  a  pointer  to 
the  waveforms  in  the  SDF.  Having  the  system  search  through  the  archived  data  in  this 
manner  rather  than  having  the  analyst  locate  and  load  the  desired  waveforms  is 
perhaps  the  principal  difference  between  the  SR1S  Pearce  subsystem  which  is  docu¬ 
mented  in  this  section  and  the  CSS  system  which  was  described  in  the  previous  section 
The  header  for  each  waveform  which  is  retrieved  from  the  SDF  is  displayed  on  the  ras¬ 
ter  screen,  and,  if  the  analyst  chooses,  the  waveform  itself  is  plotted  on  the  screen  of 
the  G-80.  The  system  then  prompts  the  analyst  either  to  add  this  new  waveform  to  the 
database  which  will  be  measured  or  to  reject  it.  If  the  waveform  is  accepted,  then  the 
system  transfers  the  waveform  to  a  memory  buffer,  and  it  retains  certain  information 
from  the  SDF  header,  such  as  the  epicenter-to-station  distance,  forward  azimuth,  and 
back  azimuth.  The  system  asks  the  operator  whether  he  wants  to  examine  more 
waveforms,  and  if  he  does,  then  it  continues  reading  the  EAF  and  retrieves  the  SDF 
record  corresponding  to  the  next  desired  (i.e. ,  P)  arrival  in  the  list.  This  process  con¬ 
tinues  until  there  are  no  more  type  1  EAF  records  for  this  event  or  until  a  maximum  of 
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30  waveforms  have  been  added  to  the  memory  buffer.  (If  more  than  30  waveforms 
exist,  the  operator  can  process  them  all  by  performing  the  analysis  on  the  first  30  and 
then  using  the  submenu  option  GETSEGWF  again  to  access  the  next  30,  and  so  on  ) 


The  GETSEGWF  option  automatically  calls  the  functions  of  the  DISPLAY  submenu 
option  to  plot  any  combination  of  up  to  18  waveforms,  chosen  from  the  memory  buffer, 
simultaneously  on  the  screen.  (It  should  be  noted,  however,  that  displaying  more  than 
8  to  10  waveforms,  depending  on  the  number  of  points  per  waveform,  usually  causes  an 
objectionable  amount  of  flicker,  and  it  greatly  reduces  the  vertical  resolution  of  each 
waveform.)  The  operator  can  subsequently  use  the  DISPLAY  option,  independently  of 
the  GETSEGWF  option,  any  time  the  IDLE  element  is  blinking.  This  will  permit  him  to 
change  the  combination  of  waveforms  which  is  shown  on  the  screen  so  that,  for  exam¬ 
ple,  he  can  examine  the  next  set  of  8  or  10  waveforms  out  of  the  maximum  of  30  which 
he  has  chosen.  DISPLAY  also  allows  the  operator  to  choose  how  many  points  of  the 
waveforms  are  to  be  plotted.  He  will  usually  select  the  entire  waveform,  since  the  P 
arrival,  and  hence  also  the  pP  and  sP  arrivals,  are  close  to  the  end  of  some  of  the  seg¬ 
mented  data  windows.  The  chosen  waveforms  are  then  plotted  on  the  G-B0.  The  traces 
are  labelled  by  a  reference  trace  number  which  will  be  used  in  further 
operator/machine  interactive  commands,  by  the  station  name,  and  by  the  source-to- 
receiver  distance  and  forward  azimuth,  Pointers  are  drawn  on  the  traces  marking  the 
arrival  times  for  the  P  phase,  which  are  calculated  using  the  Herrin  (1968)  travel-time 
tables,  the  source-to-receiver  distances  taken  from  the  SDF,  and  the  source  depth 
taken  from  the  EAF  (unfortunately,  this  value  was  always  set  to  zero  for  the  A1  events) 
The  software  for  accessing  the  travel-time  tables  and  for  computing  the  predicted  P 
arrival  times  using  them  is  the  same  as  is  used  by  the  Event  Location  System  [ELS] 
subsystem  (Teledyne  Geotech,  1982),  so  we  shall  not  describe  it  herein  In  order  to 
examine  the  waveforms,  and  particularly  the  predicted  P  arrivals,  more  carefully,  the 
analyst  may  choose  to  invoke  the  MAINMENU  option  to  return  control  to  the  IAS.  From 
the  IAS  subsystem,  he  can  use  the  main  menu  of  that  subsystem  to  expand  or 
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compress  the  plotted  trace  horizontally,  magnify  or  demagnify  it  vertically,  or  filter  it 
When  he  is  satisfied  with  the  adjustments  which  he  has  made  to  the  trace,  he  can  then 
return  to  the  Pearce  subsystem  by  again  selecting  the  SUBMENU  option  on  the  main 
menu. 

The  markers  which  have  been  drawn  on  the  waveforms  using  the  event  location 
taken  from  the  EAF  will,  of  course,  be  only  approximate  markers  to  the  true  P  arrival 
times.  The  discrepancy  between  the  predicted  and  actual  arrival  times  will  be  particu¬ 
larly  pronounced  if  the  depth  which  is  listed  in  the  EAF  is  incorrect  (as  we  have  pointed 
out.  this  is  the  case  for  the  A1  dataset,  which  lists  default  depths  of  zero  for  all  133 
events).  The  operator  is  therefore  meant  to  use  these  markers  only  as  a  guide  in 
searching  for  the  P  arrivals  visually  When  he  believes  he  has  found  a  true  signal  onset 
on  a  given  trace,  he  will  invoke  the  P-P1CKS  option  of  the  submenu,  and  then  he  will 
move  a  cursor  across  the  G-8C  screen  until  it  coincides  with  the  observed  signal  onset 
on  that  waveform  trace.  (Cursor  control  on  the  G-80  is  provided  by  a  pen  attached  to  a 
data  tablet.)  He  will  then  depress  the  pen  on  the  tablet,  transmitting  a  "hit”  to  the  sys¬ 
tem,  and  a  new  P  marker  with  the  associated  time  will  be  displayed.  The  analyst  may 
make  up  to  32  picks  to  adjust  and  readjust  the  P  phase  markers  on  all  of  the  traces 
until  he  is  satisfied  that  they  are  correct.  He  will  then  press  the  "FNI"  button  on  the 
G-80  keyboard,  and  the  predicted  P  markers  will  disappear  from  the  screen.  They  will 
be  replaced  by  new  ones  drawn  at  those  points  on  each  of  the  waveforms  which  were 
last  selected  by  the  operator  as  being  the  true  signal  start  points.  The  time 
corresponding  to  this  point  on  each  trace  is  retained  in  memory  as  the  corrected  P 
arrival  time.  The  operator  may  choose  the  P-PICKS  option  repeatedly  until  he  is 
satisified  with  the  positions  of  all  the  f’-markers.  As  a  further  guide  in  making  the 
selection  of  the  P  onset,  the  operator  may  invoke  the  T-PICK  option,  which  writes  on 
the  screen  the  time  corresponding  to  any  "hit"  made  on  the  waveform  trace  using  the 
pen  and  data  tablet  without  changing  the  value  of  the  chosen  P  arrival  time. 


The  operator  will  next  wish  to  exercise  the  DEPTH  option  on  the  submenu  This 
option  first  draws  a  horizontal  axis  on  the  bottom  of  the  screen,  marked  off  in  incre¬ 
ments  of  10  km  and  labelled  in  increments  of  50  km  from  0  to  600  km.  This  axis  will  be 
referred  to  as  the  depth  axis.  (Note  that  it  has  no  logical  connection  with  the  waveform 
traces  which  are  displayed  on  the  screen  above  it.)  The  operator  will  then  use  the  pen 
and  data  tablet  to  position  a  cursor  on  this  axis  and  "hit"  a  trial  depth.  Using  the  event 
epicenter  and  origin  time  taken  from  the  EAF  and  the  depth  which  was  chosen  on  the 
axis,  the  system  then  computes  pP-P,  sP-P,  and  PcP-P  delay  times  for  each  of  the 
stations  corresponding  to  the  displayed  waveforms.  Once  again,  the  P  travel  times  are 
computed  using  the  Herrin  (1968)  tables,  but  now  the  travel  times  for  pP  and  sP  are 
computed  using  polynomial  formulas  taken  from  the  ADAPS  system  (Teledyne  Geotech, 
1968),  and  those  for  PcP  are  computed  using  the  Jeffreys-Bullen  (1940)  tables.  The 
software  for  accessing  the  travel-time  tables  and  tables  of  polynomial  coefficients  and 
for  computing  the  travel  times  from  them  is  the  same  as  that  which  is  used  in  the  ELS 
subsystem  of  SR1S  (Teledyne  Geotech,  1982). 

Two  points  ought  to  be  noted  about  the  computation  of  these  delay  times.  First, 
mixing  travel-time  tables  in  this  manner  will  obviously  bias  the  results  since  different 
earth  models  are  being  used  to  compute  the  travel  times  of  the  different  phases  We 
also  point  out  that  the  polynomial  fit  which  is  used  for  the  sP  travel  times  is  inap¬ 
propriate  for  depths  shallower  than  about  30  km.  Second,  it  is  obviously  inaccurate  to 
choose  an  arbitrary  trial  depth  for  the  event  while  leaving  the  epicenter  fixed  to  the 
value  which  is  listed  in  the  EAF.  The  latitude,  longitude,  and  depth  variables  are  not 
independent,  so  changing  one  of  them  would  affect  the  other  two  in  the  location  algo¬ 
rithm  and  hence  in  the  computation  of  the  travel  times.  In  spite  of  these  two  objec¬ 
tions,  the  computation  of  the  delay  times  for  the  three  reflected  phases  is  adequate  for 
the  purposes  for  which  it  is  used,  as  we  shall  now  describe.  The  calculated  delay  times 
are  added  to  the  P  arrival  times  which  were  selected  using  P-PICK5,  and  markers  are 
then  drawn  on  each  waveform  trace  corresponding  to  the  predicted  pP,  sP,  and  PcP 


arrivals.  As  with  the  initial  P  markers  which  were  drawn  by  the  system  but  were  then 
replaced  by  analyst  picks,  these  markers  are  meant  simply  to  assist  the  analyst  in 
finding  the  phases  visually.  The  advantage  of  having  these  markers  is  that  they  provide 
a  means  for  relating  the  individual  traces  on  the  display  to  one  another  That  is,  the 
markers  serve  to  account  for  the  rollout  with  distance  of  the  reflected  phase  arrival 
times  in  a  manner  which  is  at  least  consistent,  if  not  strictly  accurate.  It  is  this  rollout 
which  we  wish  to  illustrate  with  these  markers.  As  Goncz  and  Barker  (1978)  have 
demonstrated,  it  is  easier  to  make  pP  picks  on  a  screenfull  of  traces  simultaneously 
than  on  each  one  of  those  traces  individually.  If  an  analyst  searches  for  pP  on  each 
trace  independently  of  the  other  traces,  he  may  end  up  picking  coda  arrivals  with 
short  delay  times,  indicating  a  shallow  source  depth,  on  some  traces  while  picking 
arrivals  with  long  delay  times,  indicating  a  deep  focus,  on  other  traces.  This  problem 
can  be  alleviated  if  it  can  be  demonstrated  that  the  pP  pick  which  the  analyst  makes 
on  one  trace  is  consistent  with  the  picks  which  he  makes  on  the  traces  for  the  other 
stations.  The  markers  which  are  drawn  by  the  DEPTH  option  are  a  means  for  illustrat¬ 
ing  this  inter-station  consistency  graphically.  The  analyst  will  use  this  feature  by 
choosing  a  trial  depth  on  the  depth  axis  and  noting  the  positions  of  the  pP.  sP,  and 
PcP  markers  which  are  drawn  on  each  trace.  If  the  depth  which  he  chose  is  the  right 
one,  a  consistent  pattern  will  appear  of  arrivals  in  the  codas  occurring  at  increasing 
delay  times  with  increasing  source-to-receiver  distance  in  the  same  manner  as  do  the 
markers  drawn  on  the  traces.  If  no  such  pattern  is  apparent,  the  operator  may  then 
"hit"  another  depth  on  the  axis  to  draw  a  new  set  of  markers.  He  will  continue  choosing 
trial  depths  either  until  he  is  satisfied  that  no  reflected  phases  can  be  detected  on  any 
of  the  traces  or  until  he  finds  a  depth  for  which  the  pattern  of  markers  seems  to  match 
certain  observed  arrivals  in  the  codas  for  at  least  some  of  the  traces.  Again,  we 
emphasize  that  this  match  will  be  only  approximate,  but  the  positions  of  the  markers 
may  be  adequate  to  guide  the  analyst  to  note  visually  the  true  onsets  of  the  reflected 
phases  In  particular,  on  those  traces  on  which  he  would  not  otherwise  note  pP,  he 
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may  now  be  able  to  detect  it  by  examining  closely  that  portion  of  the  wavetrain  sur¬ 
rounding  the  marker,  since  this  is  the  time  window  within  which  pP  would  arrive  if  the 
event  were  at  that  same  source  depth  for  which  the  predicted  pP  delay  time  does  in 
fact  agree  with  the  observed  delay  at  those  stations  for  which  he  was  able  to  identify 
the  pP  arrival.  If  he  nevertheless  is  still  unable  to  find  pfJ  on  some  traces,  the  analyst 
will  use  the  marker  to  show  him  where  within  the  coda  he  should  measure  the  coda 
amplitude  level  which  will  then  be  used  as  an  upper  bound  on  the  amplitude  of  the 
undetected  pP  phase. 

Once  the  analyst  either  has  identified  the  pP  and  sP  phases  or  has  decided,  on  the 
basis  of  the  cursor  positions,  where  within  the  P- wave  coda  he  will  make  measurements 
bounding  these  undetected  phases,  he  will  next  use  the  POLARITY  submenu  option  to 
make  picks  of  the  first  motion  direction  of  P,  pP ,  and  s P.  He  will  enter  the  reference 
number  of  the  trace  on  the  screen  for  which  he  wishes  to  make  the  picks,  and,  in 
response  to  a  prompt  appearing  on  the  screen,  he  will  enter  a  code  of  "1",  "2",  or  "3"  to 
denote  which  phase  he  is  picking.  He  then  will  enter  a  code  of  or  "0"  as  the 

polarity  pick.  The  system  will  substitute  default  values  of  "0"  for  the  polarities  of  all 
phases  which  the  analyst  does  not  explicitly  pick.  The  analyst  thus  need  pick  only 
those  phases  with  clear  first  motions  (usually  only  P  phases,  and  only  those  at  certain 
stations).  The  polarity  picks  can  be  made  in  any  order,  and  they  can  be  repeated  if  the 
analyst  recognizes  that  he  has  made  a  mistake.  As  a  guide  to  error  detection,  the 
or  "0"  code  which  the  analyst  enters  is  echoed  on  the  graphics  screen. 

The  final  set  of  measurements  is  made  using  the  submenu  option  AMPL1M.  As  he 
did  for  the  polarity  picks,  the  analyst  will  enter  information  defining  which  trace  he  is 
about  to  measure.  In  response  to  a  prompt  appearing  on  the  screen,  he  will  also  enter 
a  numeric  code  telling  the  system  whether  he  is  about  to  measure  the  amplitude 
bounds  Pnuj,  Pmta,  pPmax,  pPmin.  sPmex,  or  sP^.  These  six  measurements  may  be 
made  in  any  order,  and  they  may  be  repeated  any  number  of  times.  Default  values  of 
zero  are  supplied  by  the  system  for  any  measurements  not  made  by  the  analyst,  these 
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default  values  would  be  used,  in  particular,  for  the  minimum  amplitude  bounds  of 
undetected  phases.  The  analyst  will  then  use  the  pen  and  data  tablet  to  position  a  hor¬ 
izontal  cursor  on  the  screen  somewhat  above  (if  he  is  measuring  the  maximum  bound) 
or  below  (if  he  is  measuring  the  minimum  bound)  the  maximum  excursion  of  the  upsw¬ 
ing  peak  of  the  phase.  When  the  cursor  is  where  he  wants  it,  it  will  register  a  "hit"  by 
depressing  the  pen  on  the  data  tablet.  He  will  then  move  the  cursor  to  a  position 
somewhat  below  (for  the  maximum  bound)  or  above  (for  the  minimum  bound)  the  max¬ 
imum  excursion  of  the  downswing  peak,  and  register  a  second  "hit”.  The  system  then 
computes  the  amplitude  bound  as  being  the  difference  between  the  cursor  positions  for 
the  "hits"  on  the  upswing  and  downswing  peaks.  If  the  analyst  wishes,  he  may  instead 
measure  first  the  downswing  and  then  the  upswing  peaks,  since  the  system  uses  the 
absolute  value  of  the  difference  between  the  cursor  positions. 

In  any  event,  the  analyst  will  use  AMPL1M  to  take  the  difference  between  the  upsw¬ 
ing  and  the  downswing  "hits"  to  compute  P^,  and  he  will  use  it  to  compute  Pmln  in  the 
same  way.  Note  that  the  system  does  not  compute  the  difference  Pm&x-Pmin.  rather, 
both  of  these  values  are  written  into  the  output  file  along  with  the  maximum  and 
minimum  bounds  for  pP  and  sP  Both  the  maximum  and  minimum  bounds  are 
required,  since  the  Pearce  algorithm  uses  the  bounds  on  the  ratios  P/pP  and  P/sP 
The  minimum  bound  for  the  ratio  P/pP  is  given  by  P^/pP^,  and  the  maximum 
bound  is  given  by  Pmsx/pPm in.  with  similar  ratios  bounding  the  ratio  P/sP.  (Actually, 
in  order  to  avoid  computational  difficulties  with  large  numbers,  the  algorithm  will 
instead  use  the  reciprocal  ratios  pP/ P  and  sP/  P  if  the  original  ratios  are  too  large 
This  would  of  course  happen  when  the  P- wave  coda  is  measured  in  place  of  undetected 
surface  reflections;  in  this  case,  pPmtn  =  sPmin  =  0.)  The  f’-wave  measurements  which 
the  analyst  makes  using  the  AMPJ.IM  option  are  shown  in  Figure  15.  Similar  measure¬ 
ments  would  be  made  for  pP  and  sP,  if  they  are  detected,  or  for  the  coda  maximum  if 
these  phases  are  undetected.  The  amplitude  bounds  which  are  measured  in  this 
manner  are  reflected  on  the  raster  screen  so  that  the  analyst  can  review  them  and  can 


repeat  the  measurements  if  he  considers  them  to  be  unsatisfactory 

The  amplitude  information  is  converted  by  the  system  from  inches  on  the  C-80 
screen  to  nanometers  on  the  seismogram  using  the  scale  factor  in  the  SDK  header 
This  information  in  turn  is  taken  from  the  gain  factor  for  data  tapes  which  is  listed  in 
the  SENSOR  DAT  file  for  that  station  (c/.  DS  documentation  volume,  Teledyne  Geotech, 
1982).  It  also  uses  whatever  scaling  factors  were  applied  by  the  DISPLAY  submenu 
option  and  by  the  IAS  main  menu  features  for  waveform  trace  magnification.  It  should 
be  mentioned  that  for  the  seismograms  from  certain  stations,  the  gain  factors  in  the 
SENSOR.DAT  flies  were  in  error  by  a  factor  of  1000  since  they  were  scaled  to  microns 
rather  than  to  millimicrons;  as  we  have  pointed  out,  the  Pearce  analysis  is  unaffected 
by  calibration  problems  of  this  sort,  since  the  gain  factors  cancel  out  when  the  P/pP 
and  P/ sP  ratios  are  taken.  No  correction  for  instrument  response  is  applied  in  order 
to  convert  the  measured  amplitudes  on  the  seismogram  to  amplitudes  of  ground  dis¬ 
placement;  a  feature  for  making  this  correction  can  be  added  to  the  AMPLIM  option  at 
a  later  time  if  it  is  concluded  that  the  difference  between  the  period  of  the  P  phase  and 
that  of  thepP  and  sP  phases  is  large  enough  to  cause  a  significant  difference  between 
the  instrument  gain  factors  which  should  be  applied  to  the  measurements  of  those 
phases.  As  it  stands  now.  we  consider  the  period-related  difference  in  gain  factors  to 
be  one  of  the  many  sources  of  amplitude  variation  for  which  we  wish  to  compensate 
manually  by  placing  the  horizontal  cursors  far  enough  above  and  below  the  amplitude 
peaks  that  the  difference  Pma~Pmn  would  represent  sufficiently  large  "error  bars"  for 
the  measurement  of  the  P-wave  amplitude,  and  similarly  for  the  other  two  phases.  In 
the  next  section  of  this  report  we  shall  discuss  how  large  these  error  bars  should  be. 

It  will  be  noted  that  the  submenu  option  DEPTH  draws  markers  for  pP,  sP,  and 
PcP,  but  options  POLARITY  and  AMPLIM  both  ignore  the  PcP  phase  The  reason  that 
this  phase  is  marked  by  DEPTH  is  simply  to  assist  the  analyst  so  that  he  does  not  mis¬ 
take  this  core  reflection  for  one  of  the  surface  reflections. 
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When  the  analyst  has  finished  making  all  the  required  measurements  on  the 
waveforms  displayed  on  the  screen,  he  will  invoke  submenu  option  TAKEOFF  to  com¬ 
pute  the  takeoff  angles  from  the  event’s  hypocenter  to  the  stations  corresponding  to 
each  of  the  displayed  traces.  Input  to  this  computation  consists  of  the  epicenter-to- 
station  distances  taken  from  the  SDF  header  and  the  event  depth  corresponding  to  the 
last  "hit"  made  during  the  DEPTH  analysis.  The  computation  employs  a  fourth-order 
polynomial  fit  which  was  made  to  the  tables  of  P- wave  take-off  angles  of  Pho  and  Behe 
(1972).  These  tables  in  turn  were  derived  from  the  Herrin  (1968)  travel-time  tables 
The  coefficients  of  the  polynomial  fit  are  listed  in  Table  IV.  Linear  interpolation  was 
used  to  compute  the  take-off  angles  for  depths  intermediate  between  those  which  are 
listed  in  the  table.  As  is  shown  in  Table  IV,  we  assume  that  the  source-to-receiver  dis¬ 
tance  A  lies  within  the  range  20°  sis  100° .  We  are  unable  to  extend  this  lower  limit  on 
A  since  the  codas  of  regional  P  phases  are  contaminated  with  a  variety  of  triplicated 
arrivals  and  since  it  would  require  detailed  ray-tracing,  rather  than  just  an  average 
earth  model,  to  calculate  the  take-off  angle  of  the  P  phase  at  regional  distances.  It 
would  actually  be  better  to  use  30°  as  a  cut-off,  but  we  shall  allow  the  lower  limit  to  be 
20°  and  leave  it  to  the  discretion  of  the  analyst  whether  measurements  made  in  the 
distance  range  20°  ^A^  30°  ought  to  be  retained  or  deleted  in  each  individual  case 

Use  of  this  SRIS  subsystem  concludes  with  submenu  option  SCRSAVE,  which  writes 
into  a  file  the  station  name,  event-to-station  forward  azimuth,  takeoff  angle,  polarity, 
and  maximum  and  minimum  amplitude  bounds  for  the  data  analyzed  thus  far  This 
output  file  is  organized  so  that  each  record  contains  all  the  information  pertaining  to 
one  of  the  three  phase  pairs  P.pP,  P.sP,  or  pP.sP.  Measurements  made  for  stations 
at  distances  of  less  than  20°  or  more  than  100°  are  deleted  by  SCRSAVE,  as  are  meas¬ 
urements  in  which  the  maximum  amplitude  bound  of  some  phase  is  zero  (indicating 
that  the  phase  was  not  measured)  or  for  which  the  minimum  amplitude  bounds  of  both 
phases  are  zero  (as  is  frequently  the  case  for  the  phase  pair  pP:sP).  Default  values  of 
unity,  which  can  be  overridden  by  the  analyst  in  response  to  a  prompt  from  the 


Table  IV 


Polynomial  Fit  to  the  Tables  of  /J-wave  Take-ofT  Angles 


Takeoff  angle  =  A4  A4  +  Aa  A3  +  Az  Az  +  /t,  A  +  A0  where  : 


20°  <;  A  <  30° 

depth  (km) 

A* 

^3 

Az 

Ao 

0 

-8.9961  x  10-4 

0.08B244 

-3.1096 

45.407 

-186.74 

15 

-9.9971  x  10-4 

0.097157 

-3.3788 

48.217 

-188.12 

40 

— 1.0342 x  10-3 

0.097221 

-3.2019 

41.208 

-101.49 

100 

-1.5370x  10-3 

0.14629 

-5.0034 

70.699 

-2B7.B3 

150 

— 1.7324  x  10~3 

0.16475 

-5.6646 

61.710 

-356.31 

200 

— 1.6357  x  10~3 

0.15310 

-5.1682 

72.834 

-299.44 

250 

-1.41 15x  10-3 

0.12897 

-4.2170 

56.630 

-198.14 

300 

-1.3592X  10-4 

-1.4688  x  10"3 

0.73169 

-25  825 

311.84 

350 

-3. 6024  x  10-4 

0.020701 

-0.095277 

-11.947 

224  68 

400 

4. 2282  x  10“4 

-0.057965 

2.8319 

-59.673 

513.84 

450 

1.2112X  10~3 

-0.13608 

5.6935 

-105.52 

786.59 

500 

1.9307X  10~3 

-0.20600 

8.1988 

-144.65 

1013  3 

550 

1.9975X  10“3 

-0.20784 

8.0600 

-138.51 

952.43 

_ SQQ _ 

_ 1.4360X  IQ"3 

-0.14637 

5.5508 

-93.342 

653.70 

30®  <.  b<:  85® 

depth  (km) 

A* 

^3 

Az 

A\ 

Ao 

0 

-1.3478X  10 

3. 0859  x  10~4 

-0.025787 

0  69203 

23.8B2 

15 

-1.5862x  10-fl 

3  6270 x  10-4 

-0  030168 

0.80422 

27.251 

40 

—1.8867  x  1O“0 

4. 3052  X  10‘4 

-0.035485 

0.89992 

35  170 

100 

— 1.8528  x  10  6 

4. 2264  x  10~4 

-0.034783 

0  86392 

36.543 

150 

— 1.7764  x  1O-0 

4.0319X  10“4 

-0  032953 

0  78147 

38  643 

200 

-1.6109X  10® 

3.6568  x  10-4 

-0.029809 

0.65428 

41  627 

250 

-1.3332x  10“a 

3.0505  x  10'4 

-0.025089 

0.48568 

45.199 

300 

— 1.5685  x  10~a 

3.51 18  x  10'4 

-0.027981 

0  52881 

47  045 

350 

— 1.4794  x  10“8 

3.2904  x  10'4 

-0.025856 

0.41824 

50  795 

400 

-1.1750x  10'a 

2. 5625  x  10“4 

-0.019411 

0.14766 

56.892 

450 

-9.9141  x  10~7 

2.1119X  10~4 

-0.015206 

-0.051250 

62.402 

500 

-7.6472  x  10’7 

1.5393  x  10~4 

-9. 7699  x  10~3 

-0.30480 

68.961 

550 

-4.0815X  10-7 

6.8354  x  10~5 

-2.0155  x  10~3 

-0.64711 

76.984 

_ SCO _ 

6. 2078 x  IQ'6 

-4.8654 x  IQ-5 

6. 9045 x  IQ-3 

-1  1284 

87.341 

■7.8222 x  10  4 

0.22419 

-21.429 

697  66  i 

■8.8889  x  10~4 

0.25482 

-24  362 

793  18 

■9.9552  x  10  4 

0.28656 

-27.515 

900  94 

•1 .0134  x  10~3 

0.29098 

-27  870 

910.44 

•9.4224  x  10'4 

0.27153 

-26.106 

857  75 

■1 .0667  x  10~3 

0.30530 

-29  144 

948  92 

9.4222X  lO-4 

0.27112 

-26  029 

855.16 

•9.9558  x  10~4 

0.2B570 

-27.352 

895  79 

•1.0311  x  10~3 

0.29534 

-28  220 

922  50 

•1 .0844  x  10’3 

0.30985 

-29  532 

962  82 

■9.9599  x  10~4 

0.27567 

-26  418 

B69  50 

■9.7778  x  10-4 

0.28032 

-26  B20 

881  95 

■1.0311  x  10'3 

0.29487 

-28.139 

922  72 

■8.7112X  10~4 

0.25084 

-24  116 

601  63 

system,  are  also  written  into  this  file  for  the  pP  and  sP  reflection  coefficients  and  for 
the  ratio  of  the  seismic  velocities  at  the  source  depth  to  those  at  the  surface  (If  the 
analyst  allows  the  default  values  to  be  used,  he  will  change  them  to  more  realistic 
values  at  a  later  time  simply  by  editing  this  output  file  in  the  same  manner  as  he  would 
any  other  formatted  file  consisting  of  alphanumeric  data.)  He  may  then  request  more 
waveforms  via  submenu  option  GETSEGWF,  or  he  may  return  to  the  main  menu  of  the 
IAS  and  then  terminate  SR1S  execution 

In  the  second  (classified)  volume  of  this  report,  we  shall  describe  the  application 
of  the  SR1S  Pearce  analysis  subsystem  to  the  segmented  waveforms  of  the  Ai  dataset. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

For  explosions  with  mb  >  5,  pP  amplitude  constraints  are  sufficient  to  elim¬ 
inate  the  possibility  of  an  earthquake  at  depths  ,»  3-5  km  ,  provided  adequate 
azimuthal  coverage  exists.  The  CSS  data  set  was  inadequate  to  test  explosions 
for  mb  <  5.  Explosions  consistently  produced  P  coda  of  small  enough  amplitude 
over  a  network  of  stations,  that  they  could  not  be  mistaken  for  crustal  earth¬ 
quakes  at  depths  greater  than  3  to  5  km.  Those  earthquakes  that  failed  the 
discrimination  test  were  inadequately  recored  in  at  least  two  teleseismic  azimu¬ 
thal  quadrants.  One  station,  MAJO,  exhibited  poor  performance  with  respect  to 
other  stations  in  this  regard.  MAJO  short  period  waveforms  from  East  Kazahk 
were  characterized  by  an  emergent  P  wave  and  a  large  second  arrival  This 
behavior  is  most  likely  due  to  local  anomalous  structure  at  Matsushiro,  Japan 

Synthetic  tests  indicate  that  the  Pearce  algorithm  will  likely  eliminate  the 
earthquake  at  depth  hypothesis  if  4  to  5  stations  can  be  found  with  pP/P  ampli¬ 
tude  ratios  of  1.0  or  less,  and  the  stations  provide  270  degree  azimuthal  cover¬ 
age.  Because  of  detection  thresholds  of  the  GDSN  network,  it  was  not  possible  to 
operationally  test  whether  the  method  would  continue  to  work  at  signal-to-noise 
ratios  closer  to  1. 

Earthquakes  produced  coda  levels  and  secondary  phases  consistent  with 
double-couple  solutions  at  depths  greater  them  3  km.  Often  "constrained”  focal 
mechanisms  could  be  inferred.  Some  of  these  mechanisms  were  found  to  con¬ 
tradict  focal  mechanisms  based  on  long-period  data.  It  is  possible  that  short 
period  data  are  sensitive  to  different  and  more  local  stress  systems  than  long 
period  seismic  data.  A  concerted  effort  to  incorporate  more  polarity  data  would 
be  helpful  when  waveform  data  is  sparse  in  a  particular  quadrant 

sP  amplitudes  were  found  to  be  less  reliable  than  pP  amplitudes.  Both 


uncertainties  in  SV-to-P  surface  reflection  and  S  wave  attenuation  above  the 


source  may  be  responsible  Detailed  crustal  models  may  improve  the  SV-to-P 
reflection  modeling,  but  the  results  of  Pooley  ot  al  (19133)  indicate  that  near 
source  heterogeneity  is  also  responsible  for  deflection  of  b'V  energy  near  the 
focal  planes  Shorter  wavelength  S  waves  should  be  more  sensitive  to  any  I. acral 
velocity  gradients  than  the  longer  wavelength  P  waves  at  the  same  frequency 

Bandpass  filtering  of  digital  waveforms  was  found  to  be  useful  for  reduction 
of  noise  while  measuring  time-domain  amplitudes.  Filtering  did  not  greatly 
increase  the  detection  of  depth  phases  Both  P  coda  and  depth  phases  have 
similar  frequency  bandwidths  Consequently,  the  depth  phases  are  not  likely  to 
be  separated  from  the  coda  by  simple  filtering.  Three  component  polarization 
filtering  and  beam-forming  should  be  expected  to  improve  depth  phase  detec¬ 
tion  when  local  coda  generation  is  important.  It  was  not  possible  to  test  this 
supposition  with  this  limited  study.  RSTN  short  period  three  component  data 
should  offer  a  suitable  test  of  whether  polarization  filtering  can  improve  depth 
phase  detection  by  a  reduction  in  locally  generated  P  coda  signal 

The  Pearce  algorithm  was  found  to  be  too  slow  to  run  interactively  under 
the  present  CSS  Unix  operating  environment  Furthermore,  batch  operation  of 
the  processor  on  the  multitude  of  measurements  made  during  a  day's  analysis 
of  waveform  data  may  take  several  hours  to  complete.  This  mode  of  operation 
makes  it  difficult  for  the  operator  to  recall  all  the  detailed  logic  that  may  have 
been  involved  in  the  selection  of  specific  secondary  arrivals  from  one  day  to  the 
next.  Also,  the  seismologist  may  submit  more  than  one  alternative  set  of  ampli¬ 
tude  bounds  to  the  processor  for  later  analysis.  Access  to  a  mainframe  com¬ 
puter  operating  in  a  batch  mode  would  shorten  turn  around  time,  and  eliminate 
confusion  for  the  seismologist  attempting  to  examine  the  many  possibilities 
Some  gains  in  efficiency  in  the  Pearce  algorithm  are  possible  with  changes  in 


In  an  operational  environment,  where  the  Pearce  algorithm  may  be  used  to 
test  the  hypothesis  of  a  double-couple  mechanism  at  depth,  it  will  be  necessary 
to  quickly  and  automatically  generate  the  amplitude  bounds  from  the  seismo¬ 
grams  for  each  depth  range  of  interest,  Such  an  algorithm  should  allow  for  the 
uncertainties  in  the  initial  P  arrivals,  and  estimate  noise  prior  to  the  P  arrival 
for  comparison  with  the  P  coda.  The  seismologist  should  be  able  to  quickly  esti¬ 
mate  the  P  arrival  times  for  use  by  the  processor.  The  use  of  predicted  P  times 
based  on  a  preliminary  origin  should  also  be  available.  The  depth  interval,  and 
set  of  amplitude  bounds  could  be  quickly  stored  as  an  alternative  hypothesis  for 
the  Pearce  processor.  Such  an  automated  "pP  prospector"  would  allow  the 
seismologist  to  review  the  prospects  of  several  hypotheses  and  select  the  ’t 
likely  alternative. 
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NAME,  pc  arc  c  -  Pearce  focal  plane  analysis  program 
SYNOPSIS 

pearce  infile  [10,5]  [inconsistencies] 

DESCRimON 

Pearce  examines  a  gridded  solution  space  of  slip,  dip  and  strike  angles  for 
possible  double-couple  focal  mechanisms  consistent  with  amplitude  bounds  and 
polarities  on  P,  pP,  and  sP  phases.  The  amplitude  bounds  on  P,  pP,  and  sP  are 
converted  to  amplitude  ratio  bounds  for  either  pP/P,  sP/P,  or  sP/pP.  These 
extremal  ratio  bounds  are  considered  as  seperate  observations  and  the  bounded 
three-dimensional  solution  space  of  (slip. dip, strike)  is  searched  for  all  possible 
solutions  consistent  with  the  observations.  The  (slip, dip, strike)  space  may  be 
gridded  at  either  5  or  10  degree  intervals.  The  observations  are  provided  in  a 
formatted  input  file  xnfile.  Pearce  creates  two  files  called  printout,  and  platout. 
Printout  is  a  comprehensive  listing  of  the  processing.  Plotout  is  an  unformatted 
file  containing  results  for  input  to  fplanes,  or  histogram. 

The  second  command  line  argument  is  the  step  size  either  5  or  10  degrees, 
and  defaults  to  5  degrees.  The  third  command  line  argument  is  the  total 
number  of  inconsistenices  that  will  be  checked  for  each  respective  solution. 

This  option  can  speed  up  processing  if  the  user  is  not  interested  in  the  total 
number  inconsistent  observations  for  each  possible  solution. 

Infile  is  a  formated  input  file; 
title  (80al) 

nobservations  ppfactl  svel/pvel  (i2,2fl0.5) 
followed  by  amplitude  observations  of  the  form; 

name,  azimuth,  incidence,  Ppolarity.ampl,amp2.pPpolanty,amp3,amp4,ppfact,( _ ,sp,ra 

with  format(a8,2fl0.3,al,2fl0.3,al.2fl0.3,lf6.3,a2). 

A  sample  of  which  is  given  below  (blanks  are  specified  by  "_!'): 

10  04/06/79  18:30:05.20  41.860  77.630  10  5.10  0.00  0 
07 _ 0.7 _ 0.5 

chto _ 135.2 1 0 _ 23. 680- _ 3. 000 _ G .  000 + _ 1 4 . 000 _ 24 . 000 0  750 

chto _ 135.210 _ 23.680- _ 3.000 _ 6.000 _ 0.000 _ 24.000 _ D  500sp 

chto _ 135.210 _ 23.680 _ 14.000 _ 24.000 _ 0.000 _ 24  000 _ D.bOOra 

grfo _ 303.600 _ 2 1 . 030 _ 37. 000 _ 4 1 . 000 _ DO .  000 _ 9 1 . 000 _ 0. 500sp 

grfo - 303.600 - 21.030 _ 37.000 _ 41.000 _ 141.000 _ 175. 000 _ 0.750 _ 

majo - 75.620 _ 20.810 _ 6.000 _ 7.000 _ 13,000 _ 14.000 _ D.500sp 

majo _ 75.620 _ 20.81 0 _ 6. 000 _ 7. 000 _ 1 5. 000 _ 1 9. 000 _ 0. 750_ 

Hie  default  pP  attenuation  correction  is  ppfactl,  svel/pvel  is  the  s  wave  -to-  p 
wave  velocity  ratio  to  be  used  in  calculating  the  sP  takeoff  angle,  and  SV-to-P 
reflection  coefficient.  Amplitude  ratio  bounds  are  calculated  from  the  bounds; 
ampl/amp4,  and  amp2/amp3  with  the  polarity  information  given  by  +  _or  -. 
Attenuation  of  the  secondary  phase  may  be  indicated  by  ppfact  or  spfact.  The 
type  of  ratio  bounds,  pP/P,  sP/P,  or  sP/pP,  are  indicated  by  the  last  character 
string  on  the  line  as  either  blank  ( _ )  for  pP/P,  sp  for  sP/P,  and  ra  for  sP/pP. 

SEE  ALSO 

histagramQ  fplanesQ 


BUGS 

If  the  reduced  redundancy  check  option  is  used,  histogram  will  not  give  a 
true  picture  of  the  total  number  of  inconsistencies  for  each  observation. 
Maximum  number  of  observations  is  31,  limited  by  the  word  size  on  the  VAX. 
Pearce  is  a  shell  script  that  calls  an  f77  program,  and  the  interface  with 
UNIX  through  standard  input  should  be  changed  to  acces  the  comand  line 
argumenst  directly,  SV-to-P  reflection  coefficient  is  simple  minded  'ITie 
incidence  angle  of  SV  at  the  surface  is  assumed  equal  to  the  takeoff  angle  at 
the  source.  Also,  the  format  of  the  input  file  for  a  fortran  program  is  unfor¬ 
giving.  Watch  the  blanks  and  field  alignment. 

AUTHOR 

originally;  R.G.  Pearce  AWRE  report  No  0  4/79 

latest  modifications;  K.L.  McLaughlin  Teledyne-Geotech,  Alexandria  Va  Aug 
63 


NAMK:  fplancs  -  read  pcarcc  output  and  generate  graphical  output  of  solution 
space 

SYNOI’SIS 

fjUanes 

DNSCKHTION 

b\ planus  is  an  interactive  program  for  viewing  the  solution  space  derived 
from  the  pearce  focal  plane  algorithm.  The  user  is  prompted  for  each  input 
Three  versions  of  the  solution  space  may  viewed. 

1)  The  vector  diagram  of  f slip, dip, strike  j  space  Mach  acceptable  solution  is 
shown  as  a  vector  with  the  strike  azimuth  attached  to  a  grid  point 
representing  the  slip  and  dip. 

2)  equal  area  lower  hemisphere  projection  of  the  acceptable  focal  planes, 
with  or  without  P,T,I3  vectors.  The  stations  used  in  the  pearce  focal  plane 
analysis  are  shown  at  the  projection  of  the  P  wave  takeoff  angle  by  the  sta¬ 
tion  4  letter  code 

3)  equal  area,  or  equal  angle,  lower  hemisphere  projection  of  just  the  P.T.B 
axis  of  the  acceptable  focal  planes.  The  stations  used  in  the  analysis  are 
shown  at  the  P  wave  takeofT  angle  by  their  4  letter  code. 

The  acceptable  number  of  inconsistencies  the  user  is  willing  to  tolerate  in 
the  solution  space  is  an  option.  The  user  may  wish  to  replot  any  of  the  solution 
space  representations  with  a  change  of  the  number  of  inconsistencies  tolerated 
TTiis  option  allows  the  user  to  examine  the  robustness  of  the  solution  space  to 
loss  of  data.  If  the  pearce  analysis  was  done  with  a  grid  point  of  b  degrees,  the 
user  may  use  either  a  5  or  10  degree  grid  space  for  display  of  the  solution  space 
The  input  data  is  an  unformatted  output  file  from  histogram  and  consists  of 
a  header  followed  by  93312  integers  each  representing  a  specific  slip. dip, strike 
with  a  5  degree  grid.  Each  integer  is  the  total  number  of  inconsistencies  for  that 
possible  solution  that  were  found  by  pearce. 

SEE  ALSO 

pearce ()  histogram{) 

BUGS 

The  user  should  be  able  to  specify  which  data  elements  are  to  be  cllim- 
inated  from  the  analysis,  and  not  just  the  total  number  of  inconsistencies 
that  are  to  be  tollerated. 

AUTHOR 

K.L.  McLaughlin  Teledyne-Geotech,  Alexandria  Va  Aug  133 


NAME:  histogram  -  read  the  pearcc  output  and  generate  histogram 


SY  NOI»SIS 

histogram,  infile  outfile  [lO.b]  [yj 

DKSCRIITION 

Histogram  reads  the  plotout  file  from  pearce  and  generates  statistics  on 
the  unacceptable  solutions  for  each  possible  solution  examined  Kach  observa¬ 
tion  has  the  total  number  of  inconsistent  solutions  tabulated.  A  matrix  of 
mutual  inconsistencies  is  tabulated  to  examine  the  redundancies  with  pairs  of 
observation,  infile.  Histogram  creates  one  file  called  outfile.  Outfile  is  a  an 
unformated  file  for  input  to  fplanes  where  the  integer  array  from  pearcc  has 
been  converted  from  the  coding  of  individual  inconsistent  observations  to  the 
total  number  of  inconsistent  observations. 

The  third  command  line  argument  is  the  step  size  either  b  or  10  degrees, 
and  defaults  to  b  degrees.  The  fourth  command  line  argument  is  whether  to  do 
a  long  or  short  statistics.  Statistics  and  histogram  are  standard  output 

SKE  AISO 

pearce  ()  fplanes () 

BUGS 

If  the  reduced  redundancy  check  option  of  pearce  is  used,  histogram  will 
not  give  a  true  picture  of  the  total  number  of  inconsistencies  for  each 
observation. 

AUTHOR 

K.L.  McLaughlin  Teledyne-Geotech,  Alexandria  Va  Aug  (33 
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